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ABSTRACT 


This  study  describes  in-house  efforts  that  comprise  two  primary  objectives.  (1)  the  de¬ 
velopment  of  a  prediction  rationale  for  estimating  the  magnitude  and  character  of  gunfire 
induced,  structurally-borne  vibration  which,  in  turn,  represents  a  definition  of  equipment 
vibration  inputs,  and  (2)  the  synthesis  of  a  laboratory  vibration  teat  procedure  which 
may  be  used  to  qualify  future  equipment. 

Power  analogies,  relating  gunfire  blast  to  structural  response,  are  meshed  with  a  spatial 
parameter  to  describe  vibration  levels  that  are  functionally  related  to  the  gunfire  blast 
pulse  in  terms  of  gun  configuration,  muzzle  energy,  gunfire  rate,  and  the  vector  dis¬ 
tance  separating  the  equipment  from  the  gun  muzzles. 

A  basic,  normalized,  test  level  function  has  been  synthesized  and  integrated  into  a  labora¬ 
tory  vibration  test  procedure.  The  test  technique  developed  from  this  work  has  been  pub¬ 
lished  as  Method  T519  of  MlLrSTD-810B. 
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SECTION  I 
INTRODUCTION 


Gunfire  vibration  data  programs  and  studies  have  yielded  few,  if  any,  correlatives  with  which 
the  character  and  magnitude  of  the  resultant  structural  vibration  can  be  associated.  From  the 
equipment  viewpoint,  apart  from  protective  design  considerations,  one  desires  a  rationale  for 
estimating  vibration  levels,  spectral  shape,  and  so  on,  in  order  that  one  may  synthesize  a 
simulation  environment  to  which  equipment  may  be  exposed  and  from  which  the  successful 
emergence  will  provide  good  confidence  that  the  item  will  survive  and  properly  operate  dur¬ 
ing  real  gunfire  exposure. 

Prediction  technology,  always  important,  is  especially  so  for  the  case  of  gunfire  vibration. 
Normally,  one  may  provide  estimates  of  the  vibration  of  an  aircraft  structure  and  indulge 
the  latitude  of  assigning  a  general  and  averaged  description  over  a  zone  of  the  structure,  say 
six  feet  or  so — adjusting  the  levels  from  zone  to  zone.  In  contrast  the  large  magnitudes  and 
sharp  gradients  of  gunfire  vibration  do  not  permit  a  comparable  luxury  for  the  case  of  gun¬ 
fire  because  within  a  distance  of  six  feet  the  vibration  levels  may  change  one  order  of  mag¬ 
nitude  Moreover,  gunfire  levels,  though  ic  ■  tively  brief  in  time  duration,  are  sufficiently 
high  (as  much  as  one  or  two  orders  above  the  normal,  ambient  vibration  levels)  to  inhibit 
the  customary  practice  of  test  level  ‘''^caling-up."  Here,  we  allude  to  the  test  level  magnifica¬ 
tion  technique  employed  to  introduce  an  equivalent  laboratory  time-to-failure  —  a  technique 
which  serves  as  compensation  for  the  necessarily  limited  exposure  times  provided  by  labora¬ 
tory  tests. 

The  stress,  therefore,  is  on  a  prediction  rationale,  as  reliable  and  accurate  as  the  technology 
will  allow,  in  order  that  gunfire  vibration  levels  may  be  derived  and  real  time  simulation 
techniques  applied. 

This  work  embraces  two  elements  essential  to  the  development  and  conclusions  of  the  pro¬ 
gram.  One,  that  there  exists  a  power  relationship  between  the  gun  pulses  and  the  structural 
response  Two,  that  the  gunfire  blast  energy  is  coupled,  in  concentrated  form  near  the  gun 
muzzle  and  decays  in  some  inverse  form  of  the  distance  from  the  muzzle.  All  of  the  pro¬ 
cedures  employed,  analytics  used,  and  data  processing  techniques  introduced  reflect  either 
one  or  both  of  these  assertions. 

It  is  not  tiie  intent  of  this  introduction  to  historidze;  to  catalog  a  lengthy  list  of  equipment 
failures  and  resultant  Held  retrofit  programs — these  are  largely  a  matter  of  record  and  of 
references,  a  number  of  which  are  contained  in  this  report.  Rather,  it  is  hoped  that  the  des¬ 
cription  of  the  somewhat  unique  character  of  gunfire  vibration  and  the  identification  of  the 
major  modelling  prindples  is  suffidentiy  adumbrated  in  this  discussion  to  identify  the  basic 
dements  of  the  problem  area  and  in  so  doing,  furnish  the  guidance  necessary  for  the  de¬ 
tailed  approach  to  the  problem.  A  brief  outline  and  description  of  pertinent  details  of  this 
report  are  given  below. 

Section  II  is  concerned  with  the  energy  content  and  character  of  gun  blast  pulses  transformed 
into  the  frequency  domain  by  means  of  Fourier  Series.  Various  gun  configurations  are  ex¬ 
amined;  differences  and  commonalities  noted.  The  equations  of  the  Fourier  Series  used  in 
the  analysis  of  the  gun  blast  loading  are  contained  in  Appendix  I. 

Section  III  involves  the  apprehension,  derivation,  and  selection  of  a  broad-band,  random 
power  model.  The  relationship  between  structural  responie  anJ  pulse  spectral  power  is  shown. 
Structural  response  is  expressed  in  terms  of  gun  muzzle  energy,  tiring  rate,  and  number  of 
guns. 
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Section  IV  discusses  pertinent  details  of  gunfire  data  used  in  this  study  and  considers  the 
relationship  of  aircraft  structural  response  to  gunfire  models. 

Section  V  contains  a  regression  analysis  applied  to  four  gunfire  configurations.  Data  analysis 
is  used  to  provide  a  test  for  the  power  model.  A  description  of  the  regression  process  is  con¬ 
tained  in  Appendix  II. 

Section  VI  is  a  refinement  and  final  variation  of  Section  V  from  which  a  basic  normalized 
prediction  curve  is  derived.  An  example  is  included  to  demonstrate  the  application  of  the 
method.  A  resulting  gunfire  vibration  test  method  is  included  in  Appendix  IV. 

Section  VII  elaborates  upon  work  elements  previously  covered,  proposes  a  tentative  method 
by  which  the  gunfire  vibration  field  may  be  adjusted  as  a  function  of  the  gun  stand  off 
distance,  discusses  a  laboratory  simulation  technique  for  the  generation  and  control  of  gunfire 
low  frequency  content,  and,  finally,  explores  some  concepts  concerning  future  extension  of 
the  gunfire  model. 

Section  VIII  contains  the  conclusions. 
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SECTION  n 

PULSE  SPECTRAL  ANALYSIS 

In  aircraft  vibration  problems,  one  is  apt  to  begin  with  M.  Fourier — anticipates  ending  with 
him — and  usually  finishes  up  between  the  illustrious  Frenchman. 

To  begin  then,  a  typical  pressure-time  history  of  a  gur.blast  is  shown  in  Figure  1.  The  pulse 
duration,  t..,  is  set  at  approximately  10-J  secs — an  average  value  determined  from  micro¬ 
phone  data  taken  from  the  F-4E  and  F-6A  aircraft  [1]  [2].  The  flush  mounted  micro¬ 
phones  were  located  approximately  2-3  feet  in  front  of  the  gun  muzzles  and  in  a  zone  fre¬ 
quently  referred  to  as  the  near  field  loading  area  of  the  aircraft  structure  (Figure  8). 

Superimposed  over  the  pulse  is  its  smoothed  substitute,  the  Friedlander  blast  pulse — the  shape 
is  described  mathematically  by  the  included  expression.  The  energy  of  the  pressure-time 
excitation,  that  is  coupled  to  the  aircraft  structure  results  in  what  amounts  to  the  chief 
object  of  our  study ;  structural  vibration.  The  Friedlander  pulse  is  the  pressure-time  function 
that  we  propose  to  subject  to  Fourier  analysis.  From  this  endeavor,  it  is  hoped  that  certain 
details  of  the  spectra!  character  will  prove  sufficiently  informative  and  useful  to  assist  in 
the  construction  of  structural  response  models  and  in  the  synthesis  of  gunfire  simulation 
methods. 

Analysis  of  the  pulse  is  applied  to  three  main  cases  of  gunfire: 

Case  I.  The  Constant  Firing  Rate  Model — in  which  a  single  gun  is  assumed  to  generate  the 
Friedlander  pulse  at  a  fixed  period  and  may  be  characterized,  spectrally,  by  the  use  of  the 
Fourier  series. 

Case  II.  Variable  Firing  Rate  Model — Here,  the  gun  requires  time  to  come  up  to  constant 
speed,  stays  here  briefly,  and  then  reverses  cycles  to  the  cease  fire  state.  A  modification 
of  the  Fourier  series  is  used  as  an  approximation  of  this  model. 

Case  III.  Unsynchronized,  Multiple  Gun  Arrays — This  model  is  approximated  by  multiple 
pulse  trains,  randomly  phased  with  each  other,  and  is  characterized  by  a  special  application 
of  the  Fourier  transform. 

Each  case  is  discussed  in  the  following  paragraphs.  Mathematical  details  of  the  first  two  cases 
are  contained  in  Appendix  I.  Details  of  Case  III  are  either  discussed  in  this  Section  or  are 
referenced. 


1.  CASE  I— CONSTANT  FIRING  RATE  MODEL 

The  line  spectral  display  of  three  typical  firing  rates  (25,  50,  and  100  Hz)  is  shown  in  Figure 

2.  The  spectra  for  each  firing  rate  are  shown  with  their  respective  harmonics  separated  by  a 
frequency  interval,  Af,  equal  to  the  reciprocal  of  their  respective  gunfire  periods.  The  max¬ 
imum  value  of  the  amplitude  occurs  at  the  fundamental.  The  harmonic  amplitudes  decrease 
smoothly  with  frequency,  are  about  3  db  down  at  approximately  600  Hz,  and  from  thereon 
describe  a  slowly  oscillating  decay  rate  at  roughly  3  db/octave.  There  are  no  zero  values 
of  the  coefficient  in  this  spectrum.  One  noteworthy  property  of  the  line  spectrum,  already 
mentioned,  deserves  emphasis.  As  the  gun-fire  rate  increases,  so  in  proportion,  does  the 
magnitude  of  the  harmonics,  (Reference  1).  Similarly,  the  average  power  of  each  line  spec¬ 
trum  varies  directly  as  the  frequency  (Reference  4).  Indeed,  this  is  true  of  Fourier  series 
and  transforms  in  general.  The  power,  frequency,  amplitude  proportionality  characteristics 
is  of  some  consequence  as  will  be  seen  later  on  in  this  section  and  early  in  Secron  II. 

Taken  as  a  model  for  a  single  gun,  the  series  form  suffers  from  some  deficiencies  that  are 
worth  discussion.  What  if  the  firing  rate  fluctuates  somewhat? — and  it  does.  A  three  sigma 
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deviation  of,  say,  rtlO'y  produces  a  small  shift  in  the  spectral  lines  at  the  fundamental 
(if  =  20  Hz,  at  100  Hz)  but  the  deviation  at  the  higher  frequencies  can  be  substantia] — 
at  500  Hz  the  translation  is  100  Hz ;  at  1000  Hz,  it  is  200  Hz  and  so  on.  The  matter  is  fur¬ 
ther  complicated  if  the  firing  rate  deviation  :s  assumed  to  be  continuous  and  randomly  dis¬ 
tributed  (as  it  is,  in  fact,  usually  characterized)  because  the  Fourier  series,  in  its  applica¬ 
tion,  requires  that  the  period  be  constant. 

It  occurs  to  the  authors  that  the  firing  rate  deviation  helps  to  explain  why  accelerometers 
time  histories  of  the  aircraft  structure  (especially  noticeable  above  approximately  400  Hz) 
exhibit  amplitude  distributions  that  are  predominantly  random  in  character  rather  than 
sinusoidal.  Too,  it  is  suspected  that  this  is  reinforced  by  the  tuning-detuning  effect  on  air¬ 
craft  structures  arising  from  fluctuating  flight  loads. 

In  any  case,  the  model  has  decided  limitations — for  the  low  frequency  portion  of  the  spec¬ 
trum  it  seems  the  more  germane  and  useful.  For  the  high  frequency  region,  the  application 
appears  restricted.  Not  mentioned  thus  far  is  the  formidable  simulation  difficulty  associated 
with  the  generation  and  control  of  some  20  to  80  spectral  magnitudes. 

2.  CASE  II— VARIABLE  FIRING  RATE  MODEL 

The  variable  firing  rate  model  is  synthesized  to  reflect  the  response  time  characteristics  of 
guns  that  require  a  finite  time  interval  to  come  up  to  normal  firing  rate.  Guns  of  the 
rotating  barrel  class  (Vulcan)  and  guns  that  have  been  exposed  to  low  temperature  (Refer¬ 
ence  5)  are  typical  examples  of  the  variable  firing  rate  model.  Oscillogram  data  (during  gun¬ 
fire  initiation)  from  a  gun  of  the  Vulcan  class  indicated  a  sharp  rise  in  the  rate  to  approx¬ 
imately  50  Hz,  followed  by  a  more  gradual  increase.  This  rise  continued  for  approximately 
7/10  second  at  the  end  of  which  the  normal  firing  rate  of  100  Hz  was  reached.  From  the 
firing  history,  a  simplified  rate-time  function  was  constructed.  For  constructional  simplicity 
and  computational  convenience,  the  7/10-second  time  interval  was  rounded  off  to  1  second. 
A  total  cycling  history  consists  of  a  cycling  up  (Part  1),  a  constant  rate  (Part  2),  and  a 
cycling  down  (Part  3),  each  of  1-second  duration.  The  history  is  shown  in  Figure  3.  Seg¬ 
ments  of  the  corresponding  pulse-time  history  are  shown  in  Figure  4. 

A  full  sequence  of  bursts  is  treated  as  the  new  function,  p(t).  Periodicity  is  introduced  by 
allowing  the  burst  of  pulses  to  be  repeated  every  3  seconds.  This  amounts  to  3  second  gunfire 
bursts,  repeated  every  three  seconds. 

Coefficients  of  the  Fourier  serin  for  the  pulse  train  were  obtained  from  a  special  arrange¬ 
ment  of  the  inverse  Fourier  transform,  the  details  of  which  are  contained  in  Appendix  I. 

The  line  spectrum  results  are  shown  in  Figure  5  except  that  frequency  amplitudes  are  in¬ 
serted  every  100  Hz ;  all  others  are  omitted  for  shape  comparison  with  Case  I.  The  frequency 
content  in  between  is  necessarily  omitted  for  clarity  of  presentation.  Recall,  that  in  order  to 
present  the  assembly  of  pulses  in  terms  of  the  Fourier  series  the  entire  pulse  assembly  be¬ 
comes  p(t)  and  this  is  made  periodic  by  constructing  a  sequence  of  bursts,  periodically 
spaced.  The  choice  of  three  seconds  between  pulse  bursts  is  arbitrary — any  reasonable  inter¬ 
val  could  have  been  chosen.  Since  3  seconds  was  selected,  the  spectral  lines  are  separated  by 
approximately  1/3  Hz — a  very  dense  line  spectrum ;  hence  the  reason  for  the  spectral  omis¬ 
sion  in  Figure  5.  Probably,  the  most  interesting  feature  of  Figure  5  arises  from  the  obser¬ 
vation  that  the  cycling  process  has  not  appreciably  changed  the  main  outline  of  the  line 
spectrum.  Differences  occur  for  the  most  part  in  the  region  below  800  Hz  and  these  varia¬ 
tions  do  not  seem  noteworthy. 

Characteristic  distributions  of  side  band  frequencies  about  the  main  spectral  magnitudes 
are  shown  in  Figure  6.  Again,  for  purposes  of  clarity,  spectral  lines  are  placed  at  spaced 
intervals  to  allow  the  insertion  of  a  sample  of  the  side  bands  that  are,  in  turn,  spaced  ap- 
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proximately  1/3  Hz  apart.  Of  course,  similar  side  band  distributions  occur  about  the  omitted 
harmonics  (200,  300,  500,  600,  700  Hz  and  so  on). 

Other  gunfire  oscillograms  of  the  same  type  gun  Bhow  roll-up  times  shorter  than  that  of  the 
gun  selected  for  this  study.  Apparently,  our  gun  (gas  operated)  is  of  a  sluggish  species.  This 
note,  apart  from  being  a  commentary  on  the  variability  of  guns  (their  firing  and  cycling 
rates)  provokes  the  thought  that  if  a  shorter  roll-up  time  is  introduced  into  the  cycling  model 
then  one  would  anticipate  that  the  weighting,  provided  by  a  proportionately  longer  time  spent 
at  constant  rate,  would  result  in  a  spectrum  that  tends  to  converge  to  the  100  Hz  model  of 
Case  I. 

In  summary,  the  cycling  model  shows  no  significant  change  in  the  envelope  of  the  line  spec¬ 
tra,  which  is  the  Fourier  transform  P(f)  for  a  continuous  spectrum,  and  results  in  a  genera¬ 
tion  of  spectral  side  bands  distributed  about  the  fundamental  firing  frequency  and  about  the 
harmonics. 

A  study  of  accelerometer  time  histories  for  this  type  of  gun  shows  the  same  upper  and  low 
frequency  anomalies  noted  in  Case  I.  That  is,  amplitude  distributions  above  approximately 
400  Hz  show  a  predominant  random  character.  Whereas,  further  down  the  spectrum,  har¬ 
monic  content  appears  increasingly  dominant. 

3.  CASE  III— UNSYNCHRONIZED,  MULTIPLE  GUN  ARRAYS 

Multiple,  randomly  phased,  pul9e  trains  introduce  the  transform  aspect  of  M.  Fourier.  But, 
before  proceeding  to  this  alternative,  we  include  an  additional  assumption  comprising  this 
model  that  is  extracted  from  Case  I.  We  require  that  the  pulse  train  of  each  gun  be  approx¬ 
imately  periodic  but  allowed  to  deviate  about  its  respective  mean  firing  rate  with  a  distribu¬ 
tion  that  is  approximately  Gaussian.  From  Case  II,  if  appb'cable,  the  roll-up  and  roll-down 
effect  on  the  pulse  sequence  may  also  be  included  in  the  train  history.  Finally,  the  pulse  train 
duration,  T,  is  assumed  to  be  much  larger  than  the  time  duration  of  the  pulse,  t.. 

The  available  model  found  most  realistic  in  simulation  and  tractable  in  application  comes  from 
electrical  theory. 

The  following  discussion  outlines  the  main  elements  of  a  random  pulse  train  model.  For  more 
detailed  information,  see  Reference  4. 

Consider  a  train  of  gunfire  pulses  occurring  randomly  at  times  t„  t»  t> . . .,  t.  (Figure  7a). 
The  pulses  may  overlap  or  may  be  widely  spaced. 

The  random  pulse  train  may  be  characterized  in  the  frequency  domain  by  application  of  the 
real  translation  theorm  to  the  Fourier  transform  P,  («),  of  the  pulse  train,  as  follows: 

P,  (.)  -  P(„)  [e“  imt'  +  e“  j-ti  +  . . .  +  e"  j“ts ]  (1) 

where  P(«)  is  the  transform  of  the  pulse,  p(t).  And  the  energy  density  spectrum  E (.)  is 

E(«)  -  ^*|P(-)|J[e  ~imil  +  e  “  iuti  +  ...  +  e  “  j-ts]*  (2) 

where  t„  t-,  t„  t*  are  pulse  occurrence  times  and  are  distributed  randomly  over  an  ensemble 
of  N  pulses. 

The  exponential  series  of  the  right  side  of  Eqation  2  is  treated  as  a  two-dimensional  random 
walk  problem  in  which  Equation  2  emerges  in  terms  of  the  power  density  spectrum,  f  (») . 

*<-)  «r*P*(o)*(.)  +  r|P(«) |J«l/2»  (3) 

where  8  («)  is  the  transform  of  the  unit  impulse  function,  8  (t) 

and  8(«)  *  1  when- *  0 
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Figure  7.  Random  Pulse  Trains  (a)  and  Power  Density  Spectrum  (b) 


12 


t  («*)  "0  when  m  *  0 

The  term,  i^P^O),  represents  the  DC  power  spectrum  at «  “  0  and  is  zero  elsewhere.  Fig¬ 
ure  76  depicts  f  (•)  for  the  pulse  trains. 

When  the  DC  term  is  dropped,  Equation  3  simplifies  to 

*(•)  «r|P(.)p.A.  (4) 

If  the  energy  density  function,  |P(«)  jJ,  of  the  Friedlander  pulse  is  replaced  with  the  corres¬ 
ponding  |P(«)  I1  for  a  sawtooth  pulse,  then  Equation  4  provides  the  model  for  the  voltage 
power  density  output  (voltsJ/Hz)  of  a  random  noise  generator — a  device  that  utilizes  a 
noisy  diode  (shot  noise)  for  the  spectrum  source. 

Strictly  speaking,  the  model  described  by  Equation  4  is  not  a  precise  match  to  the  model 
requirement  of  Case  III.  The  chief  difference  resides  in  the  concept  of  random  pulse  train*. 
The  gunfire  pulse  trains  are  randomly  phased  with  each  other,  but  though  variable,  each 
train  does  exhibit  its  o  ,n  periodicity.  The  model  described  by  Equation  4  admits  of  no  pe¬ 
riodicity  whatever — only  a  statistical,  average  rate.  Differences  in  the  spectral  details  -of. 
these  two  models  have  not  been  explored  analytically  but  are  expected  to  be  the  subject  of 
future  inquiiy.  Experimental  evidence  does  exist  (Reference  6),  however,  and  the  results 
(for  two  guns  firing)  show  accelerometer  responses,  noticeably  random  in  character  begin¬ 
ning  at  the  high  frequencies  and  extending  down  to  below  200  Hz.  Even  the  response  am¬ 
plitudes  at  the  fundamental  firing  frequency  (a  beat  modulation)  can  be  more  accurately 
characterized  as  random-like  rather  than  sinusoidal. 

The  random  noise  model  then  seems  an  acceptable  choice  for  Case  III.  In  retrospect  we  note 
that  all  three  cases  possess  a  significant  and  highly  useful  commonality.  Above  approximately 
400  Hz,  all  exhibit  structural  response  behavior  that  may  be  classified  as  predominately 
random-like. 

It  is  this  feature  that  we  exploit  throughout  the  greater  part  of  the  study,  beginning  with 
Section  III. 
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SECTION  III 


THE  RESPONSE  MODEL 


To  construct  a  mode!  relating  the  muzzle  blast  pressure  to  the  resultant  structural  response, 
we  introduce  a  gun  configuration  with  the  muzzle  situated  a  perpendicular  distance  h  from 
the  aircraft  surface  (Figure  8).  The  maximum  induced  structural  response,  due  to  coupling 
with  the  blast  pulses,  is  conceived  to  be  concentrated  at  the  near  field  loading  area ;  and  as 
we  sight  down  the  distance  vector  D,  we  visualize  the  response  magnitudes  to  propagate 
radially-decaying  as  some  function  of  D.  Here,  a  restriction  is  imposed  on  the  allowable 
magnitude  of  h,  thus,  constraining  the  parameter  as  a  variable.  From  a  look  at  Figure  9, 
we  conclude  that  the  pressure  isobars  decrease  appreciably  with  respect  to  the  muzzle  near 
field  values,  when  h  is  allowed  to  be  greater  than  25  calibers.  For  our  impending  model 
then,  it  is  required  that  h  be  less  than  approximately  25  calibers.  Additional  restrictions 
come  to  mind  while  we  are  on  the  subject  of  model  assumptions,  and  it  is  here  appropriate  to 
list  others  thought  to  be  pertinent  to  the  synthesis  of  the  conceptual  model. 

1.  The  response  contribution  of  the  shock  wave  traversing  the  structural  surface  is  pre¬ 
sumed  to  be  imperceptible,  beyond  the  near  field. 


2.  The  influence  on  the  response  due  to  details  of  the  reflected  wave  fronts  and  their  inter¬ 
action  with  the  incident  wave  are  considered  to  have  been  obscured  by  the  averaging  proc¬ 
esses  of  the  response  model. 

3.  The  gun  muzzles  are  essentially  clean ;  that  is,  feature  no  diffusion,  deflectors,  or  special 
blast  attenuation  shrouds. 

4.  For  the  multiple  gun  case,  the  spacing  between  muzzles  is  restricted  to  approximately 
30  calibers. 

First,  one  must  get  from  the  gun  muzzle  to  the  structure  surface.  To  do  this,  we  consider  the 
relationship  between  the  power  spectrum  of  the  blast  pulse  and  the  power  spectrum  of  the 
response. 

We  make  use  of  energy  models  by  invoking  Parsevals  theorem  which  states  that  the  energy, 
E,  contained  in  a  nonperiodic  pulse,  is  equal  to  1/2  »  times  the  area  described  by  the  square 
of  the  absolute  value  of  the  pulse  energy  density  spectrum  |P(H)|2,  or 

E-  V4*  f|P(«)Pd-  (5) 


The  pulse  energy  per  unit  frequency  is  designated  E(«),  so 

E(«y-'%»|P(-)i*  (6) 

The  average  power  per  gun,  II,  generated  by  a  succession  of  pulses  is  equal  to  the  number 
of  pulses,  N,  divided  by  the  time,  T,  consumed  by  the  pulse  bursts,  or 


but 

N 

r,  the  average  gunfire  rate  (pulses/sec) 
so 

n  »  r,«E 


(7) 
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And,  in  spectral  power  terms 

n  (•)  -  r,«E  (•)  (8) 

Multiplying  the  energy  density  function  ifo  |P(-)|J  by  r„  we  obtain  the  spectral  power, 
n(«),  in  terms  of  |P(«)|2.  Thus 

n  (•)  -r,»  |P(«)|2  •  y%*  (9) 

We  note  that  in  Equation  9,  if  we  set  r  ■  y  where  t  is  the  time  period  average  (of  one 

pulse)  of  the  energy  density  spectrum,  the  right-hand  term  then  is  the  power  density 
spectrum,  f(«). 

Substitute  in  Equation  9 


#U)-r,.|P(.)|*.  (10) 

and  using  Equations  6,  9  and  10,  we  have 

E(„).r,-*(.)  (11) 


and  the  total  power /gun  of  the  pulse  is 


n  “  E  •  r,  - 


For  n  guns. 


n  m  E  •  r,  •  n 


n  /,< 


+  • 

f  («)  d*» 


(12) 


(13) 


To  briefly  recapitulate,  we  have  expressed  the  pulse  blast  of  the  gunfire  in  terms  of  its 
Fourier  transform  in  the  power  mode.  The  primary  reason  for  introducing  the  power  form 
is  that  we  now  propose  to  linearilly  relate  the  acceleration  power  spectral  density  of  the  air¬ 
craft  structural  response  to  the  pressure  power  spectral  density  generated  by  the  gunfire 
pulse  trains. 


If  we  knew  bow  much  of  f  («)  is  coupled  to  the  structure  and  if  one  could  describe  the 
transfer  function  |H(«)  |J  relating  power  input  to  structural  response,  G(-),  then  we  could 
designate  the  power  conversion  efficiency  as  $  (-)  and  describe  file  structural  response  thus: 


J *(-)  H 


at)  |H(«)  |2<L»  OC 


f G(»)da 
*  —  «• 


Unfortunately,  we  cannot  readily  assign  functional  values  to  *(•),  although,  as  will  be 
seen  in  Section  IV,  we  do  have  a  somewhat  better  (though  stfll  generalized)  notion  of  the 
character  of  |H(«)  |2. 


Our  approach  to  this  difficulty  is  to  make  use  of  similitudes,  equate  coupled  power  of  the 
pulses  to  structural  response,  and  eliminate  intractable  properties  as  shown  in  the  following 
steps.  Let  0(«)  |H(«)|3  f  («)  equal  the  spectral  power  flow  into  the  structure  and  let  the 
resultant  power  spectral  density  of  the  structural  response  be  a— lgm*h>il  G(«).  Equating 
spectral  power  input  to  structural  response. 


*(.)  |H(-)|2*(.)  -/«(•) 


(14) 


Substitute  Equation  11  in  14, 


*(„)  |H (•)  |J  E  («*)  rj  ™  pG(m) 

Assume  a  reference  gun  (with  zero  subscripts) , 


(15) 
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(1«) 


M-)  |H„(«)  I1  E„  (.)  •  r.  -  pG.(m) 

Dividing  Equation  15  by  16  and  rearranging  gives 

G(-)  -  E(m)/E..(»)  •  r./r.  •  #’(.)/«.(•)  *  G.<«)  •  |H(-)|*/H.  <-)|a  <17> 

The  power  transfer  efficiency,  fl(-),  and  the  power  transfer  function,  |H(«)|J,  are  assumed 
to  be  the  same  for  both  guns.  Thus,  0  («)  ”  M»)  and  |H(.)  ”  |H,(«)  |.  Expression  17  then 
becomes 

G  (-)  -  E  (-)  /E0  (-)  •  r,/r,  •  G,  (-)  (18) 

In  terms  of  overall  mean  squared  response.  Equation  18,  is  shown  as 


j G(»)d«  “  E/E.  •  r,/r.  •  jo*  (<• 


We  introduce  the  muzzle  energy  E'  by  proportioning  it  to  E 
E'  ■  «E  and  E«'  ”  a„E, 

The  proportionality  constant,  «,  is  assumed  to  be  the  same  for  both  guns 
(*  "  «•)  • 

Substituting  Equation  20  into  19  gives 
!•+•  /*+• 


^ G(«)  dm  **  E'/E'.  •  r,/r,  G.(«)  d«* 


Finally,  since  the  responses  are  normalized  to  one  gun,  we  insert  n  and  a,  to  account  for  mul¬ 
tiple  gun  configurations,  see  Equation  13.  Thus  Equation  21  becomes 


Jct(m)  dm  -  E'/E' 


•  n/n.  •  r,/r.  I  G.(»)  d» 


J  *  n 
The  appearance  of  the  gun  muzzle  energy  is  welcome.  Normally  this  property  is  readily  avail¬ 
able  and  so  provides  us  with  a  much  needed  parameter  to  facilitate  the  prediction  process. 
The  muzzle  energy  is  equal  to  mvJ/2;  where  m  is  the  projectile  mass  and  v  is  the  muzzle 
velocity. 

We  have  expressed  the  predicted  overaD  mean  squared  response  of  an  unknown  »un  conftgm- 

ation  in  terms  of  its  muzzle  energy,  the  number  of  guns,  and  the  firing  rate— all  normalized 
to  some  known  gun  system  having  like  variables. 

To  establish  the  necessary  relationship  of  the  parameters  of  Equation  22  to  their  Proper 
Aw  in  the  regression  analysis  beginning  with  Section  V,  we  now  turn  to  Section  “ 
this  section  we  discuss  general  characteristics  of  aircraft  structurnl  response,  details  of  data 
processing  and  presentation,  gunfire  configurations,  ballistic  parameters,  selection  p*  *cccfer- 
ometers.  and  review  other  pertinent  details  that,  in  their  varied  ways,  are  associated  with 
tiie  development,  evaluation,  and  application  of  the  gunfire  power  model. 
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SECTION  IV 

DATA  PROCESSING;  DISPLAY 

L  RESPONSE  SPECTRUM 

As  suggested  earlier,  the  spectral  form  of  the  response  has  influenced  the  manner  in  which 
we  have  selected,  organized,  and  presented  the  data  for  regression  analysis — an  examina¬ 
tion  of  the  response  characteristics  is  necessary  to  explain  how  the  data  has  been  processed 
and  why  it  appears  in  the  form  that  it  does,  and  in  what  way  the  response  envelope  enters 
into  the  prediction  model. 

If,  for  example,  we  select  as  an  input  Fourier  spectrum  F(f),  obtained  from  the  blast  pulse 
(Figure  76)  and  supply  a  generalized  frequency  response  function  H(f),  typical  of  aircraft 
structures  then  for  the  response  function,  R(f) : 

R(f)  «F(f)H(f) 

Phase  information  is  assumed  to  be  random;  a  continuous  random  spectrum  is  also  as¬ 
sumed;  further  F(f)  and  R  (f)  is  expressed  in  1/3  octave,  mean  squared  density  ideom: 

[H(f)P-[P(f)]*[H(f)]J  (FigurelO) 

The  resultant  spectral  response  is  shown  in  Fig  10c.  [R(f)P  exhibits  the  familiar  roll-up, 
peak,  roll-down  form  so  frequently  encountered  during  1/3  octave  data  analysis  of  acceler¬ 
ometer  outputs  in  response  to  gunfire  excitation. 

To  continue,  the  spectral  shape  of  the  response  is  generalized  by  flattopping  the  peak  region 
and  (using  log-log  axes)  representing  the  roll-up  and  roll-off  portions  of  the  curve  with 
straight  lines.  The  curve  will  then  appear  as  shown  in  Figure  11. 
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Figure  10.  Generalized  Vibration  Field  of  Aircraft  Structure 
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Figure  11.  Generalized  Response  of  Aircraft  Structure 

An  examination  of  the  gunfire  data  in  this  study  (see  Figures  41,  42,  43,  and  44,  Appendix 
III)  results  in  the  estimation  of  corner  frequencies  and  slopes  as  indicated  in  Figure  12. 
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Figure  12.  Corner  Frequencies  and  Roll-Off  Slopes  of  Structural  Response 


Now,  if  we  concern  ourselves  primarily  with  the  1/3  octave  response  between  315  to  1250 
Hz,  then  for  a  given  accelerometer  located  a  distance  D,  from  the  gun  muizle,  the  response 
data  might  appear  as  shown  in  Figure  13. 
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Figure  13. 1/3  Octaves  of  Response 

Finally,  if  we  rotate  the  plot  90°  such  that  the  frequency  axis  now  runs  into  the  page  and 
further,  if  we  include  the  responses  from  other  accelerometer  samples  located  at  their  respec¬ 
tive  d™**™*—  from  the  gun,  then  we  might  have  the  following  display  of  the  data  points. 
Figure  14. 
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Figure  14.  Presentation  of  Response  vs  Distance 


The  response,  [R(f)p,  is  displayed  as  a  function  of  the  distance  parameter,  and  the  solid 
line  represents  the  regression  line  generated  from  a  least  squares  fit  for  the  data.  Now  sub¬ 
stitute  G,  (f)  in  place  of  [R(f)]1.  This  display  of  the  data  will  be  encountered  exclusively 
in  the  following  sections.  Note  that,  for  the  regression  technique  employed  here,  the  structural 
response,  G,(f),  is  confined  to  the  plateaued  segment  of  the  spectral  response  and  that  the  re¬ 
sponse  data  is  presented  in  1/3-octave  acceleration  power  spectral  density  form. 


We  return  to  Equation  22,  which  is  restated 


/ 


+  • 

G(f)df 


E'/E',  •  n/n*  •  r,/r. 


/ 


+  - 

G.(f)df 


(22) 


We  assume  the  generalized  response  function  (Figure  10)  to  be  invariant  in  shape  but 
having  a  maximum  ordinate  value,  G„(f),  that  is  a  linear  function  of  power,  thus: 


(f)  df  -  jBG,(f)  and 


/' 


+  - 

G(f)df 1 


pG(f) 


Substituting  the  above  equations  into  Equation  22  and  rearranging  into  normalized  form, 
we  have 


G(f)  -E'/E'.«n/.«r,/r..G.<f) 


(23) 


We  sort  out  the  reference  part  of  the  power  model : 

E',  •  r.  •  n.  -  0G,(f)  (24) 

From  hereon  our  efforts  will  involve  a  two-pronged  course:  to  evaluate  and  define  the  para¬ 
meters  of  Equation  24  and  to  attempt  the  description  of  the  decay  function  relating  G„(f) 
and  the  distance  parameter,  D.  Both  efforts  involve  regression  analysis  applied  to  gunfire 
data  obtained  from  four  different  gun  configurations  of  approximately  equal  muzzle  energies. 


2.  GUN  CONFIGURATIONS 

The  response  data  from  accelerometers  located  in  two  different  aircraft  was  processed.  The 
data  was  obtained  from  three  configurations  of  the  F-4  and  one  of  the  F-fiA. 

The  nose  configuration  of  the  F-4E  is  shown  in  Figure  15,  and  is  referred  to  in  the  data  and 
charts  as  simply  F-4.  This  arrangement  features  the  M-61  type  cannon  arranged  as  six 
barrels  firing  consecutively  during  revolution  and  is  sometimes  identified  as  the  Vulcan  or 
the  gattling  gun. 

The  remaining  two  gun  configurations  consist  of  20-mm  cannons  located  in  pod  packages 
that  are  mounted  beneath  the  aircraft  surface  and  on  centerline.  They  are  designated  the 
SUU-16  and  the  MK-IV  pods  and  are  shown  in  Figure  16.  The  MK-IV  pod  consists  of  two 
barrels  aide  by  side  firing  simultaneously,  while  the  SUU-16  pod  features  one  barrel  only. 
Both  of  these  configurations  hereafter  will  be  referred  to  as  the  MK-IV  and  SUU-16. 

Two  M-S9A2  cannons  are  fuselage  mounted  in  the  top  side  of  the  F-5A  aircraft  and  forward 
of  the  cockpit  as  shown  in  Figure  17.  The  guns  are  spaced  approximately  24  inches  apart  This 
configuration  will  be  referred  to  as  the  F-5.  One  important  feature  of  the  F-5A  must  be 
stressed:  there  are  two  retractable  deflectors  located  at  the  gun  muzzles.  These  devices  are 
employed  during  gunfire,  therefore  during  the  data  acquisition  phase,  and  serve  to  deflect  the 
gaseous  charge  up  and  away  from  the  inlets  of  the  engine  intakes  in  order  to  prevent  engine 
stall.  This  sjiedal  feature  and  its  influence  on  the  decisions  concerning  the  final  stages  of 
regression  analysis  will  reappear  in  Section  VI. 

In  summary,  all  configurations  represent  variations  of  the  20-mm  cannon  and  all  feature  the 


23 


26 


Figure  17.  F-5A  Gunfire  Configuration  and  Accelerometer  Locations 


same  muzzle  energies.  They  differ  chiefly  in  the  number  of  guns  and  firing  rates.  The  two 
parameters  (number  of  guns  and  firing  rate)  figure  importantly  in  the  gunfire  model,  will 
be  referred  to  frequently,  and  are  listed  in  table  I  for  reference. 


TABLE  I 

GUNFIRE  MODEL  PARAMETERS 


Configuration 

Firing  Rate 
r(Hz) 

Number  of 
Guns,  n 

n  •  r 

(ft -lbs) 

F-4 

90 

1 

90 

39,600 

8UU-16 

100 

1 

100 

39,600 

MK-IV 

85 

2 

70 

39,600 

F-5 

27 

2 

54 

39,600 

3.  ACCELEROMETER  LOCATION 

Wherever  possible,  accelerometer  data  were  eeketed  from  transducer  locations  most  repre¬ 
sentative  of  the  primary  structure  and  positioned  as  dose  to  the  skin  surface  as  possible.  In 
genera],  the  volume  enclosed  by  the  skin  surface  to  a  depth  of  10  to  12  inches  included  mo3t  of 
the  accelerometer  locations. 

4.  DETERMINING  D 

The  D  parameter  represents  the  vector  distance  from  the  gun  muzzle  point  to  the  selected 
accelerometer  location  and  was  calculated  from  the  three  orthogonal  distances  referenced  from 
the  butt  line,  fuselage,  and  the  water  line  coordinates.  Figure  18  represents  a  typical  spatial 
arrangement  from  which  D  is  determined. 

5.  VIBRATION  DATA  CONVERSION 

All  of  the  gunfire  data  used  in  this  study  was  obtained  from  References  1,  2,  7,  and  8  and 
appeared  in  %  octave  form.  Data  not  processed  originally  in  RMS  g  terminology  was  so  con¬ 
verted.  The  g  values  were  squared  and  divided  by  their  appropriate  V4  octave  bandwidths  and 
presented  as  G„(f)  in  units  of  acceleration  power  spectral  density. 

6.  GUNFIRE  FLIGHT  CONDITIONS 

Both  F-4  and  F-6  gunfire  data  were  recorded  during  similar  flight  condition  at  approxi¬ 
mately  5000  feet  altitude  and  at  Mach  0.85. 

7.  REGRESSION  ANALYSIS 

Regression  analysis  involved  the  method  of  the  least  squares,  the  details  of  which  arc  con¬ 
tained  in  Appendix  I.  The  regression  lines  of  the  figures  are  bounded  by  95%  confidence 
limits  for  the  line  as  a  whole.  Stated  another  way :  if  each  of  the  gunfire  runs  were  repeated, 
there  exists  a  95%  probability  that  a  regression  line,  obtained  from  the  new  sample,  will 
fall  within  the  limits  indicated  by  the  confidence  band. 
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SECTION  V 

POWER  MODEL  EVALUATION 

Figures  19  and  20  show  an  example  of  the  response  amplitude  plotted  against  D  and  repre¬ 
sents  a  gunfire  vibration  sample  presented  in  both  the  linear  and  logarithmic  plane,  respec¬ 
tively.  Both  data  plots  illustrate  the  nature  of  the  decay  that  suggests  an  exponential 
model  and  so  the  basic  expression  for  the  decay  model  is  assumed  to  be  a  regression  line  of  tlie 
form: 

G.(f) -/Ur1  (26) 

For  a  given  muzzle  energy,  the  response,  G„  ( f )  is  also  proportional  to  the  firing  rate,  r,  and 
the  number  of  guns,  n  —  if  we  recall  the  power  model  derived  earlier  and  represented  by 
Equation  24.  Thus,  from  a  power  model  viewpoint.  Equation  25  can  also  be  expressed  in  the 
normalized  mode  as  follows : 

G.(f)/nn,  -  /JJD’Vmi*  (26) 

Ideally,  if  the  power  model  reigns,  we  could  normalize  the  data  of  each  configuration  (F-4, 
F-6,  SUU-16,  MK-IV)  to  their  respective  n's  and  r'a  and  since  they  have  the  same  muzzle 
energies  then  the  prediction  model  represented  by  Equation  29  requires  that  the  corresponding 
G*(f)’s  of  each  configuration  be  equal.  In  fact  (assuming  that  each  configuration  curve  has 
the  same  decay,  or  slope  value,  the  curves  would  converge  to  one  line.  But  in  the  real  work), 
(as  the  lamentation  goes)  where  at  best,  models  represent  good  approximations  and  where 
data  scatters  remains  perhaps  the  only  certainty,  we  are  prudent  to  expect  less  dramatic 
results. 

Suppose  we  compare  the  data  in  the  following  way :  pool  the  data  of  the  four  separate  config¬ 
urations,  Figure  21,  and  then  conduct  a  regression  analysis  of  the  total  assembly  of  data. 
The  regression  line  obtained  from  the  least  squares  fits  Figure  22  is,  for  this  case,  presumed 
to  be  independent  of  the  power  parameters  n  and  r  and  is  of  the  form  of  Equation  25.  Next, 
the  data  of  each  configuration  is  normalized  to  its  respective  n’s  and  r*s  (obtained  from  Table 
I)  then  the  data  is  pooled,  and  following  this,  subjected  to  a  regression  analysis  (Figure  23). 
In  comparison  to  the  unnormalized  cue,  we  would  expect  (if  the  power  model  represents  a 
viable  premise)  that,  on  the  average,  the  data  of  each  normalized  configuration  would  tend 
to  converge.  If  so,  then  the  goodness  of  fit  or  multiple  correlation  coefficient,  Cr,  would 
increase  somewhat,  the  dispersion  and  the  standard  error  S/  should  likewise  decrease.  Com¬ 
parative  values  of  C.  and  S,'  resulting  from  this  procedure  are  listed  in  Table  II. 


TABLE  II 

REGRESSION  ANALYSIS 


c. 

3/ 

Slope  (z) 

Pooled 

.6162 

.7062 

•2.64 

Pooled  and  Normalized 

.6572 

.6881 

-2.87 

As  noted  in  Table  II,  C,  and  S/  increase  and  decrease  respectively,  in  agreement  with  the  prior 
hypothesis.  Note  that  when  stated  in  the  converse,  the  hypothesis  is  not  necessarily  true.  To 
rephrase:  the  stipulation  that  Cr  increase  and  S/  decrease  is  necessary  but  not  sufficient  to 
validate  the  power  model.  Conceivably,  one  could  obtain  similar  results  using  some  other 
assumed  model.  Nevertheless,  the  improvements  in  C,  and  S,',  though  small,  are  in  the  right 
direction.  For  additional  checks  of  the  power  model,  we  must  look  elsewhere. 
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During1  the  gunfire  program  of  the  F-B  (Reference  2),  Vs  octave  vibration  data  was  re¬ 
corded  from  an  accelerometer  with  both  guns  firing  and  then  with  only  one  gun  operating. 
Here,  is  a  case  for  the  power  model  in  which  n  changes  from  2  to  n  *  1.  Under  these  con¬ 
ditions,  the  power  form  for  the  single  gun  requires  that  G„(f)  be  halved.  Figure  24  shows 
a  response  plot  for  both  cases.  On  the  whole,  the  vibration  level  does  decrease,  and,  on  the 
average,  a  decrease  of  3db  would  be  a  reasonable  estimate.  Unfortunately,  a  similar  rec¬ 
ord  was  not  available  permitting  a  corresponding  check  of  the  firing  rate,  r. 

In  summary,  the  use  of  comparative  analysis  utilizing  regression  techniques  supported  by 
separate  evidence,  in  which  the  test  variable  is  n,  lends  support  to  the  power  model. 

We  proceed  to  the  next  step  which  involves  further  refinement  of  the  data  and  the  de¬ 
termination  of  the  decay  slope  of  the  composite  (pooled)  data. 
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Figan  24.  Structural  Response  for  One  and  Two  Guns 
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SECTION  VI 

SPECIFICATION  SYNTHESIS  AND  APPLICATION 


1.  SYNTHESIS 

Earlier,  in  Section  III,  it  had  been  remarked  that  the  F-5A  employed  special  gun  blast  de¬ 
flectors  during  normal  gunfire  operations.  And  now  that  we  examine  Figure  25  and  ob¬ 
serve  that  the  slope  of  the  F-5  data  is  in  sharp  contrast  to  the  others,  (a  slope  of  approx¬ 
imately  one)  we  wonder  if  there  exists  an  association.  An  examination  of  several  accel¬ 
erometer  gunfire  records  (there  were  only  several  available)  recorded  first  with  the  de¬ 
flectors  retracted  and  then  extended  show  that,  on  the  average,  the  vibration  levels  were 
reduced  when  the  deflectors  were  deployed.  These  accelerometer  locations  were  near  to 
and  in  front  of  the  guns  and  correspond  to  most  of  the  F-5  data  points  in  the  region  cov¬ 
ered  by  D  less  than  40  inches.  Were  these  levels  in  Figure  25  to  be  increased,  the  end  ef¬ 
fect  would  be  to  increase  the  slope  of  the  regression  line. 

Whether  or  not  this  explanation  fully  accounts  for  the  slope  difference  cannot,  with  such 
limited  data,  be  determined  but  it  certainly  appears  to  be  a  contributant.  Another  factor 
may  involve  the  separation  distance  between  the  guns.  At  24  inches,  the  guns  are  about 
30  calibers  apart  and  the  near  field  loading  area,  relative  to  that  assumed  for  one  gun, 
is  now  being  distributed  over  a  more  extensive  region  of  the  aircraft  structure.  With  a 
broader  distribution  area,  one  would  expect  a  somewhat  more  gradual  decay  of  the  re¬ 
sponse  levels. 

Initially,  one  is  reluctant  to  remove  the  F-5  data;  for  one  thing,  we  will  have  reduced  the 
total  sample  size  and  keyed  the  remaining  data  to  a  single  aircraft  structure  and  to  one 
data  acquisition  and  reduction  process.  But  then  again,  being  reduced  to  one  structural 
sample  is  not  without  some  advantage.  The  coupling  and  response  data  may  be  amena¬ 
ble,  in  the  future,  to  normalization  techniques.  That  is,  we  may  be  able  to  apply  surface 
(or  volume)  density  weightings  to  other  aircraft  structural  configurations  to  obtain  an 
improved  estimate  of  G(f) — especially  for  those  structures  that  markedly  deviate  from 
what  might  now  become  the  F-4  norm. 

On  balance,  then,  and  primarily  because  of  the  use  of  deflectors  which  is  in  clear  viola¬ 
tion  of  model  assumptions  (with  consequences  that  seem  to  be  reflected  in  the  data)  we 
are  persuaded  to  set  aside  the  F-5  data  sample. 

If  the  F-5  data  is  removed,  how  will  C,  and  S/  behave  if  we  repeat  the  former  compara¬ 
tive  regression  process?  The  results  from  this  inquiry  are  shown  in  Table  III,  and  a  data 
display  of  each  case  is  shown  in  Figures  26  and  27. 


TABLE  m 

REGRESSION  ANALYSIS  WITHOUT  F-5  DATA 


c. 

s. 

Slope(z) 

Pooled,  normalized  without  F-5  data 

.7143 

.6806 

-4.58 

Pooled,  without  F-5  data 

.7190 

.6675 

-4.55 

Similar  to  earlier  results,  there  has  been  some  improvement  in  the  multiple  correlation  coeffi¬ 
cient  and  in  the  sample  standard  error  when  the  data  model  is  viewed  in  the  power  (normal¬ 
ized)  mode.  From  here  on,  we  will  utilize  the  pooled,  normalized  for  (with  the  F-5  data  re¬ 
moved)  as  the  nucleus  for  developing  a  basic  normalized  reference  level,  G»(f)  /r„n„. 
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Figure  25. 
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The  next  step  involves  a  vertical  translation  of  the  regression  line  by  approximately  2  S,'  to 
obtain  93 coverage  of  the  data  points.  But  before  this,  we  wish  to  adjust  the  slope  some¬ 
what.  Earlier,  it  had  been  observed  that  the  slope  of  the  MK-IV  data  sample  exhibits  a 
noticeably  steeper  value  than  that  of  the  F-4  and  SUU-16  samples  (see  Figure  21).  Since  we 
wish  to  be  conservative  and  also  since  the  vibration  levels  of  the  MK-IV  accelerometer  at 
D  *  78  inches  appear  to  be  sufficiently  low  to  allow  a  suspicion  that  this  data  channel  might 
have  been  either  faulty  or  unrepresentative,  we  rotate  the  regression  line  of  Figure  27  counter¬ 
clockwise  to  the  full  range  allowed  by  the  upper  and  lower  confidence  limits  of  the  line  as 
a  whole  (Figure  28). 

The  rotated  regression  line,  with  a  slope  of  approximately  —3,  is  then  translated  vertically 
to  provide  an  approximate  2  S/  coverage  of  data  points.  Further  conservatism  is  obtained 
during  rotation  by  not  rotating  the  line  about  D  mean,  as  is  customary,  but  by  moving  the 
center  of  rotation  to  a  D  of  approximately  38  inches.  This  is  done  to  insure  that  the  line 
gives  better  coverage  of  those  data  points  located  close  to  the  gun  muzzle — the  region  where 
most  gunfire  failures  are  apt  to  occur. 

The  last  step  consists  of  fairing  the  curve  into  an  asymptote  to  approximate  the  near  field 
roll-off  characteristic  and,  finally,  a  reverse  curve  is  introduced  at  the  far  end  of  D  to  ap¬ 
proximate  the  gradual  blending  of  the  gunfire  levels  with  ambient  vibration  field.  Figure  29. 

G„(f)  tends  to  be  somewhat  lower  at  D  behind  the  gun  muzzle  when  compared  to  G„(f)  at 
an  equivalent  D  in  front  of  the  muzzle  (a  regional  differentiation  indicated  in  Figure  29  by 
(-)  and  (  +  )  signs,  respectively).  The  difference,  approximately  3  db,  is  derived  from 
somewhat  scanty  data  extracted  from  the  F-6  and  F-100  samples  (References  9  and  13) 
and  so  the  minus  curve  represents  a  rather  coarse  estimate.  However,  the  existence  of  a 
minus  curve  does  have  a  rational  basis;  the  maximum  blast  pressure  field  is,  of  course,  lo¬ 
cated  directly  in  front  of  the  muzzle. 

The  final  curve  (Figure  29)  defines  G„(f)  /r-n,,  as  a  function  of  D  and  will  be  used  hence¬ 
forth  as  the  normalized  reference  level  by  which  an  estimated  level,  G(f),  is  established 
when  E,  r,  n,  and  D  are  known  or  can  be  estimated. 

2.  APPLICATION 

The  following  example  is  provided  to  demonstrate  the  application  of  the  power  model  when 
it  is  desired  to  estimate  the  vibration  levels  of  an  aircraft  equipment  exposed  to  gunfire. 

An  equipment  is  located  forward  of  two  20-mm  cannons  (n  ■  2).  The  distance  parameter 
is  determined  to  be  55  inches  (D  *  55  inches).  The  muzzle  energy  is  known  to  be  43,500 
ft/lbs  (E'  =  43.5  x  103).  The  firing  rate  is  25  Hz  (r  ■  25).  The  basic  normalized  curve 
of  Figure  29,  or  Figure  519-1  of  Appendix  IV,  is  derived  from  guns  with  an  E'.  =  39.6  x 
103  ft/lbs. 

Steps: 

1.  Refer  to  Figure  29  or  519-1.  Select  from  the  positive  curve,  G(f)/r.n,  “  10-*,  at  D  =*  55 
inches. 

2.  Substitute  the  appropriate  values  into  Equation  23. 

G(f)  -  E7E’.»nr«G(f)/rA 

-ww-(2)  ,25)  (IIH) 

-  (1.1)  (50)  (10-*) 

G(f)  ■  0.55  g*/Hz 

3.  Refer  to  Figure  519-4  of  Appendix  IV.  Assign  to  GMO  the  value,  0.55  g*/Hz. 
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Figure  28.  Rotation  and  Translation  of 
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In  Section  I,  we  noted  the  desirability  of  utilizing  real  time  as  the  ciioice  for  the  gunfir*  test 
duration.  During  this  study  information  was  collected  concerning  the  number  of  rounds  fired, 
number  of  missions  per  total  rounds  expended,  mission  duration,  gun  firing  rates,  and  so  on. 
Reports  and  data  were  collected  from  SEA  activities,  gunnery  range  sources  (Nellis  AFB), 
and  from  private  contractors.  Although  reports  were  frequently  incomplete,  an  estimate  of 
30  minutes  gunfire  exposure  time  was  determined,  on  the  average,  to  be  equivalent  to  at  least 
300  missions.  Thus,  a  period  of  30  minutes  was  chosen  as  the  basic  test  time  for  T519  of 
M1L-STD-810B  and  appears  in  the  time  schedule  of  Procedure  II  of  that  method. 

For  additional  examples  and  details,  see  Method  TS19,  MIL-STD-810B,  which  is  contained 
in  Appendix  IV. 
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SECTION  VII 

FURTHER  CONSIDERATIONS  CONCERNING  GUNFIRE 

The  following  resume  is  included  to  further  elaborate  on  matters  covered  in  the  main  body 
of  the  report,  to  discuss  simulation  methods  which  involve  the  low  frequency  portion  of  the 
gunfire  spectrum ;  and  to  outline  study  areas  the  investigation  of  which  should  prove  fruit¬ 
ful  for  the  future. 

1.  MUZZLE  DEFLECTORS:  BAFFLES 

In  earlier  assumptions,  deflectors  or  special  baffles  were  disallowed.  In  practice,  most  gun 
configurations  feature  some  form  of  blast  deflection.  The  SUU-16,  MK-IV,  and  F-4E  are 
no  exceptions  to  the  rule.  A  small  lip  or  protuberance  over  the  muzzle  and  a  metallic  webbing 
between  muzzles  are  features  of  the  SUU-16  and  MK-IV,  respectively.  A  perforated  metal 
housing  appears  around  the  F-4E  configuration.  Such  arrangements  do  not  appear  to  be 
particularly  effective,  judging  from  the  results  of  other  investigators  (Reference  1)  in 
consequence,  these  nominal  features  are  treated  here  as  if  the  gun  muzzles  are  clean — such  is 
not  the  case  with  the  F-5  configuration  as  noted  in  previous  sections  of  this  work. 

2.  MUZZLE  ENERGY  AND  BLAST  ENERGY 

In  the  early  sections  of  this  work  we  set  the  energy  of  the  blast  proportional  to  the  muzzle 
energy  (Equation  20).  This  was  done  to  establish  the  blast  power  which  was  then  equated 
to  the  power  spectrum  of  the  structural  response.  The  muzzle  energy  is  also  a  convenient 
property  consisting  of  the  muzzle  velocity  and  projectile  mass— quantities  readily  obtained 
for  a  large  variety  of  weapon  classes. 

Recent  studies  (Reference  10)  suggest  another  energy  definition  of  the  gun  blast  energy 
that  under  certain  circumstances  may  represent  an  improvement  This  model  equates  the 

blast  energy  (Ep)  to  the  difference  between  the  energy  of  the  explosive  charge  (Er)  and 

(mv2 1 

— g- 1.  Stated  in  symbols: 


where 

Ep  =  blast  energy 
Er  ”  explosive  charge  energy 
m  *  project .e  mass 
v  ~  projectile  velocity 

There  are  difficulties  with  this  model  centering  chiefly  around  different  energy  values  of 
various  explosives  (values  depending  upon  whether  pressure  or  temperature  increases  are 
used  to  determine  the  charge  energy)  and  the  relative  scarcity  of  published  information  of 
blast  data  that  is  also  related  to  the  quantity  or  type  of  propellants  used.  Nonetheless,  the 
model  merits  further  investigation — revealing  a  particular  advantage  for  the  case  when  gun 
barrel  lengths  are  varied  substantially.  For  example,  there  are  cases  in  which  the  cannon 
barrel  is  reduced  to  one-third  of  original  length  with  a  consequent  reduction  of  muzzle  veloc¬ 
ity  and  an  increase  of  blast  pressure.  Note  that,  for  this  case,  the  difference  model  predicts  an 
increase  in  blast  energy — the  muzzle  energy  model,  in  contrast,  requires  a  decrease.  In  future 
efforts  we  expect  to  study  this  model  and,  if  appropriate  values  of  Er  can  be  accumulated, 
consider  seriously  its  adoption. 
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3.  MUZZLE  DISTANCE 


The  SUU-16  and  MK-IV  pods  are  mounted  underneath  the  aircraft  surface.  The  vertical 
distance  (h)  from  the  muzzle  to  the  skin  is  approximately  12  inches,  or  15  calibers. 
Both  units  are  also  attached  to  the  aircraft  by  means  of  a  centerline  adapter.  In  this  case,  h 
is  increased  to  22  and  26  inches,  respectively.  In  calibers  (c)  this  now  represents  distances 
ranging  from  approximately  23  to  33,  values  somewhat  larger  than  .our  earlier  restriction 
of  25  calibers.  The  vibration  data  used  in  this  program  represents  a  mix  of  data  taken  with 
and  without  the  adapter.  Information  on  blast  pressure  field  variations  as  a  function  of 
distance  from  a  structural  surface  is  not  extensive,  but  Westine  (Reference  10)  con¬ 
ducted  a  series  of  measurements  which  may  be  of  considerable  value.  A  weapon  is  placed 
over  a  flat  sheet  containing  six  flush  mounted  microphones  arranged  in  a  line  perpen¬ 
dicular  to  the  trajectory  (Figure  SO) .  The  weapon  is  translated  parallel  to  the  surface,  thus 
allowing  measurements  to  be  made  at  different  horizontal  distances  (L)  from  the  gun 
muzzle  to  the  microphone  array.  Also  upon  each  positioning  of  L,  the  height  (h)  was 
varied;  for  example  at  an  L/c  of  20,  pressures  were  recorded  for  h/c  ratios  of  16,  27, 
38,  and  50. 

Reflected  pressure,  taken  an  a  function  of  the  perpendicular  distance  (h)  between  the 
muzzle  to  the  plate  surface,  was  extracted  from  the  data  of  Westine.  A  spatial  pressure 
average  was  obtained  over  a  rectangular  surface  area  (see  dashed  rectangle  of  Figure 

30)  located  forward  of  the  gun  muzzle,  beginning  at  L/c  9  20,  extending  forward  to  L/c 
9  50 ;  and  bounded  on  each  side  of  centerline,  at  ±  16  inches.  This  region  represents  an 
approximation  to  the  near  field  blast  region  where  blast  coupling  is  presumed  to  gen¬ 
erate  the  most  severe  vibration  response.  The  average  pressure  of  the  surface  was  norm¬ 
alized  to  its  maximum  value  and  is  presented  as  the  dependent  variable  of  h/c  (Figure 

31) .  The  dashed  line  has  been  added  to  the  figure  in  order  to  extrapolate  the  function  into 
a  greater  and  more  useful  region  of  h.  Finally,  the  extended  line  is  faired  into  an  inverse 
square  law  curve  beginning  at,  approximately,  h/c  >  70. 

From  this  curve,  two  conclusions  result:  (1)  The  earlier  height  restriction  imposed  on 
the  model  (h/c  <  25)  seems  justified.  (2)  The  vibration  data  of  this  model  should  be 
adjusted  to  compensate  for  the  mean  blast  pressure  drop  somewhere  near  the  region  of 
h/c  >  40.  Since  the  data  from  the  SUU-16  and  the  MK-IV  consists  of  a  composite  of  h/c 
~  22  and  h/c  ~  32,  major  adjustments  of  G„(f)  are  not  anticipated.  Nevertheless,  it 
should  be  interesting  to  review  this  aspect  and  to  determine  if,  as  seems  likely,  correlar 
tion  improves  upon  data  adjustment.  Of  course,  the  end  results  arrived  at  here  should 
be  tempered  with  qualification  that  there  exists  a  relative  paucity  of  h/c  data  and  that, 
granted  these  results,  it  remains  true  that  aerodynamic  flow  and  aircraft  surface  curva¬ 
ture  and  irregularities  are,  presumably,  not  reflected  in  the  pressure  vs  h/c  data  used 
here.  Additional  studies  involving  in-flight  data  acquisition,  in  which  h/c  is  made  a  var¬ 
iable,  is  to  be  much  desired. 

In  the  meantime,  we  propose  to  use  the  curve  of  Figure  32  as  in  interim  guideline.  Ac¬ 
cordingly,  the  curve  has  been  squared  and  transformed  into  a  power  model  which  allows 
for  a  reduction  of  the  acceleration  power  spectral  density  of  the  structural  response 
G„  (f),  or  G„,.„  in  the  terminology  of  Method  T519.  This  curve  will  be  appended  to  the 
standard  in  the  immediate  future. 


4.  LOW  FREQUENCY  CONSIDERATIONS 

In  Section  IT,  we  noted  that,  for  the  case  of  a  single  cannon  firing  at  an  approximate  fixe* 
rate,  the  character  of  the  vibration  response  seems  most  realistically  simulated  by  in 
corporating  a  fundamental  frequency  component  including  harmonics  that  extend  to 
cross-over  region  at  approximately  400  Hz.  For  the  case  of  multiple,  unsynchronized  gum 
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Gan  fiut  vi  Distance  Parameter,  h/c 


low-frequency  random  characteristics  appear  predominate  and  this  is  reinforced  as  the 
num!>er  of  guns  increase,  as  in  the  case  of  the  four  gun  configuration  of  the  F-100  (Ref¬ 
erence  9).  In  the  first  case,  it  is  tempting  to  consider  an  array  of  pulsed  decaying  sin¬ 
usoids  api>ended  to  the  low  frequency  segment  of  the  random  spectrum.  For  the  second 
case,  low  frequency  extension  of  the  basic  random  spectrum  down  to  the  fundamental  fir¬ 
ing  frequency  ap|>ears  acceptable.  Clearly,  the  first  case  poses  the  greater  difficulty  and 
so  will  be  discussed  further. 

It  should  be  possible  to  disable  the  desired  low  frequency  section  of  most  random  shaker 
systems  (by-pass,  roll-off,  remove  equalizer  cards  or  preamplifiers,  etc.)  and  to  insert  the 
required  number  of  spectral  components  up  to  approximately  400  or  500  Hz.  Deviations 
about  the  mean  firing  rate  and,  possibly,  simulation  of  the  gun  runup-rundown  varia¬ 
tions  might  be  achieved  through  frequency  modulation  of  a  ramp  or  pulse  generator,  as 
shown  in  Figure  33.  Admittedly,  unforseen  difficulties  may  be  encountered  with  such  an 
approach;  therefore,  follow-up  laboratory  studies  are  planned  to  more  thoroughly  ex¬ 
plore  this  and  other  possible  techniques. 


5.  GUNFIRE  MODEL  EXTENSIONS 

It  is  of  interest  to  consider  some  extensional  possibilities  of  the  basic  gunfire  model.  To 
briefly  review,  the  present  model  results  from  an  attempt  to  provide  a  tractable  descrip¬ 
tion  of  the  vibration  field  of  the  outer  structure  including  the  contiguous  primary  struc¬ 
ture  of  the  aircraft.  Thus,  accelerometer  data  of  the  interior  locations  (usually  proximate 
to  or  on  black  box  or  other  equipments)  was,  with  purpose,  withheld.  The  D  vector  to  be 
consistent  with  this  model  must  then  terminate  at  the  structural  surface,  or  more  pre¬ 
cisely,  in  the  region  of  the  surface  bounded  by  a  depth  of  approximately  10  inches.  In  prac¬ 
tice  and  in  the  absence  of  extensional  notions  of  the  interior  vibration  field,  we  presently 
extrapolate  the  vector  to  the  desired  point.  Perhaps  this  is  sufficient — to  put  it  plaintively, 
the  best  that  can  be  done.  But  it  has  been  observed  by  the  authors  that  as  the  vibration 
propagates  into  the  interior  structure,  the  spectrum  appears  to  undergo  transformation 
that  (if  these  changes  can  be  shown  to  possess  at  least  a  statistical  consistency)  may  be  ex¬ 
ploited  to  further  useful  purpose.  That  is,  in  a  general  way,  the  high  frequencies  are  attenu¬ 
ated  and,  the  low  frequencies  relatively  speaking,  undergo  emphasis.  In  brief,  there  is  a  spec¬ 
tral  downshift  of  the  maximum  responses;  seemingly  at  the  expense  of  the  high  frequency 
amplitudes.  This  seems  to  occur  because  the  interior  dynamics  of  equipment  bays  are  charac¬ 
terized  by  concentrations  of  high  density  equipments  with  tower  stiffness  to  mass  ratios  than 
are  usually  found  near  the  aircraft  surface.  In  transit  damping  tosses  and  inertial  loading,  no 
doubt,  add  their  contributions. 

Given  such  a  hopeful,  viable  view  of  vibration  propagation  (coupled  with  a  dynamicists  last 
resort,  statistical  acceptability)  we  may  then  allow  the  D  vector  to  terminate  at  or  near 
the  structural  surface  and  add  an  R  vector  normal  to  the  surface,  extending  into  the  equip¬ 
ment  bay  to  the  desired  location.  Figure  34  illustrates  the  notion. 

The  earlier  view  of  G«(f)  is  of  the  general  functional  form; 

G,(f)  «f(nrEp/D) 

and,  in  the  extensional  view,  might  appear  as 
G»(f)  -f(nrE,/DR) 

The  significant  parameter  changes  (envelope)  of  the  altered  spectrum  might  be  keyed 
functionally  to  R. 
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Combined  Sinusoidal  and  Random  Simulation  of  Gunfire 


SECTION  VIII 
CONCLUSIONS 


Structural  vibration,  induced  by  gunfire,  is  shown  to  be  characterized  by  maximum  levels 
concentrated  in  front  of  and  near  the  gun  muzzle.  The  levels  exhibit  a  sharply  decaying  be¬ 
havior  as  a  function  of  the  distance  from  the  muzzle  and  can  be  approximated  by  an  ex¬ 
ponential  function. 

It  is  shown  that  a  power  model  relating  the  gun  parameters  (firing  rate,  muzzle  energy, 
number  of  guns)  to  the  vibration  levels  provides  a  viable  tool  for  the  estimation  of  vibration 
levels  and  facilitates  the  derivation  of  gunfire  test  techniques. 

A  vibration  test  method  has  been  synthesized  that  utilizes  a  basic  random  vibration  spectrum. 
For  the  case  of  single  gun  configurations,  it  is  desirable,  in  the  future,  to  determine  a  suitable 
simulation  method  by  which  the  low  frequency  region  (below  300  Hz)  may  be  augmented 
by  vibration  content  consisting  of  a  frequency  modulated  fundamental  and  its  harmonic 
series. 
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APPENDIX  I 


HARMONIC  ANALYSIS  OF  GUN  BLAST 
PULSES  ON  AIRCRAFT  STRUCTURES 

In  the  analysis  of  aircraft  structural  vibration  excited  by  the  aerodynamic  loading  associated 
with  gunfire,  it  is  necessary  to  obtain  the  relative  magnitude  of  the  various  harmonic  com¬ 
ponents  of  this  forcing  function.  This  is  required  when  solving  for  the  transfer  function 
giving  the  relation  between  the  aerodynamic  power  input  and  the  resultant  vibratory  power 
in  the  aircraft  structure.  This  aerodynamic  loading  is  a  result  of  the  shock  wave  and  the 
associated  overpressure  region  impinging  against  the  outer  surface  of  the  aircraft  structure 
during  its  propagation  from  the  gun  muzzle.  The  variation  in  pressure  as  a  function  of  time 
at  any  point  on  the  aircraft  structure  subjected  to  this  gun  blast  loading  may  be  obtained 
from  the  following  equation  from  (Reference  10) : 


p'(t>  -p*<l-£>  *e-i 


1.  CONSTANT  FIRING  RATE  MODEL 

A  representative  pattern  of  this  gun  blast  loading  at  a  constant  firing  rate  is  shown  in 
Figure  35.  This  is  the  basic  configuration,  or  mathematical  model,  of  the  gunfire  pattern. 
This  constant  firing  rate  may  be  expressed  by  a  periodic  function,  p(t),  which  can  be 
represented  by  the  following  Fourier  series  in  the  complex  form. 

The  Fourier  series  in  complex  form  for  the  periodic  function,  p(t),  shown  in  Figure  36  is 


+- 


P(t)  -  2 


k  -  — • 


A  +  T 


where  a*  *  — 


dt 


The  ak  is  the  complex  Fourier  coefficient,  T  is  the  period  of  f(t),  and  t»  is  arbitrary.  It  is 
convenient  to  choose  t.  =  0  so  that  the  limits  of  integration  are  0  and  T. 

pi(t)  is  defined  over  the  interval  0  <  t  <  T  as  follows: 


P  (t)  *  Pw  •  ^1  -  •  e  t.  o<t<t. 


P  (t)  -  0 


t,  <  t  <T 
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The  expression  for  a,  which  must  be  evaluated  according!)’  is  : 


dt 


This  integral  can  be  evaluated  analytically  and  the  final  expression  for  ak  will  be  put  in  the 
polar  form  of  a  complex  number. 


Let  «*  —  2* 


k 

T 


^  Jo  Tt< t/ » 


Let  Pk 15 +  j«* 


Since  J* e**dx  *4  e“  and  ^ xe“  dx  =  p-  (ax  -  1) 


(Pm  e-*1  _  Pm  e~»* 
Tb.  ~  Pk  Tt.  p*i 


'-"‘‘-"It 


a*  «  (<1  +  P*t)  e-r*‘-p,t.e-^)  £ 

-1+pA) 


Since  p*  =  r-  +  M 
% 


|  ('-,l  +  1+ *•«-  « + «■»  ) 


*>  ■»-£  (ny  ’  (*•*• + j-t.) 


•* " (r+ktr)*  *  (^e**'t’ +  608  “  i 8ln  (•**•)) 
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Let  (1  +  j.j*t0)  *=  Vl  +  e1"*,  where  ok  13  arc  tan  (-kt„) 

t  e*2-4  /  V 

K  =  Pm  *  ^  [1+  (mX)*]  (C0S  +  J  ^e0,kt”  “  S'n 

The  complex  number  enclosed  in  brackets  will  be  put  in  polar  form.  Denote  this  number 
by  hk 

hk  -  cos  («kt„)  +  j  (e<uxt0  -  sin  (^t.) 
hk  =  VcosJ  (.h,t„)  +  [e„,kt0  -  sfnT»i£TP  e*8‘ 


hk  =  VT+le«kt„)J -2(e»»t0)  sin  (»kt0)  eJ»‘ 

.  _  _  .  /e«kt„  —  sin  («kt..)\ 

where  ft -arc  tan  (  C08  ) 

Substituting  this  expression  for  hk  in  the  last  expression  for  ak  gives 

_  t„  V 1  +  <e«kt„)J-2  (e«ktn)  sin  UO  _  J  (ft  -  2ak) 
a‘"P''eT - rr+UU’i - - 


Let  Ak  =  p„, 


t„  V  1  +  (e«kt„)J  -2  (e^U  sin  («kU 
ET  [1  +  Ut,)*] 


Let  *k  53  ft  -  2«k 

Then  the  final  polar  form  for  the  complex  Fourier  coefficient  is 


a*  —  Ake 


j*k 


Substituting  this  expression  for  ak  in  the  series  gives 


P(t) 


i*kej-kt 


The  Ak  is  the  amplitude  of  the  harmonics  at  negative  as  well  as  positive  frequencies.  The 
amplitudes  for  positive  frequencies  only  can  be  obtained  by  changing  from  this  complex  form 
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to  the  corresponding  cosine  form  of  the  series  which  is: 


p(t)  =  A„  +  y^2Ahcos  (»tt  +  4>k) 

•  k  - ;  '  ;  y 

This  change  in  form  can  be  accomplish^  as  follows: 

p(t)  =  ^Akej  <",kt  +  =  A0  +  ^  (AktP  *  *k)  A.k  e  i  <•**  +  ♦»> 

k-  «  k - t 

Since  A  =  Ak,  <f>  =  -  4>k  and  «  =  -  »k 

-4c  — k  — k 

p(t)  -  A.  +  ^ Ak  (ei (wkt  +  *>.)  +  e  -j M  +  *„) 

k  -1 


P(t) 


A0  + 


^2  Akcos  (»kt+  *k) 

k  -1 


A,„  the  DC  term,  is  the  average  value  of  p(t). 

The  amplitude  of  the  harmonic  at  the  frequency  c*  =  2*  k  radians 
per  unit  of  time  is  2Ak.  T 


The  complex  Fourier  coefficient,  ak,  from  the  preceding  equations,  has  been  computed  by  the 
use  of  a  computer  program,  for  constant  firing-rates  of  25,  50,  and  100  rounds/sec.  The  mag¬ 
nitude  is  plotted  in  Figure  2  as  a  function  of  the  harmonic  frequencies. 


2.  VARIABLE  FIRING  RATE  MODEL 

A  schematic  figure  of  the  pressure  pulses  for  a  variable  firing  rate  mathematical  model  is 
shown  in  Figure  4.  The  period  T,  of  the  firing  rate  is  plotted  versus  the  occurence  time,  t, 
as  shown  in  Figure  3.  The  variation  of  the  instantaneous  blast  pressure,  p,  as  a  function  of 
time  is  determined  from,  the  following  mathematical  analysis: 

Let  the  sequence  of  N  pressure  pulses  shown  in  Figure  36  be  represented  by  the  non-periodic 
function  p(t)  which  is  defined  as  follows: 

p(t)  *  6,  t  <  t, 

a _ t.  \  -<•  —  «»> 

P(t>  =p„,  <i_S5z±2)«— —  V.<t<  (tk  +  t.) 

for  k  “  1, 2, 3, - ,  N 

0,  (tk  +  t.)  <t<tk  +  I,k  - 1,2,3,—,  (N—l) 

0,  t  >  (tN  +  t.) 


P(t)  = 
P(t)  - 
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The  duration  of  each  pulse  is  t..  and  tk  is  the  occurence  times  of  the  leading  edge  of  the  pulses. 
Thc  Fourier  transform  of  p(t)  us 


g  (-)*=  A<-)e^“> 

AM  and  <p{~)  are  the  amplitude  and  phase  spectruma  of  p(t),  respectively. 
Since  p(t)  **  0,  for  t<  t,  and  t>  ts  +  tf 


This  integral  can  be  expressed  as  the  sum  of  N  integrals  since  p(t)  =  0  between  adjacent 
pulses. 

*<->  -  £  (  jtl  (1-^)  e=^«-i-*dt 

Changing  the  variable  of  integration  in  each  integral  by  the  substitution  z  *r  t  — t*  in  the 
kth  terms,  g (»)  becomes : 
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This  result  shows  that  the  transform  of  p(t)  is  the  product  of  the  transform  of  a  single 
pulse  occurring  at  zero  time  and  a  factor  which  consists  of  a  sum  of  N  complex  numbers 
each  of  unit  magnitude  and  with  a  phase  angle  which  is  determined  by  one  of  the  occurrence 
times,  t*. 


Since  e 


cos  («tii)  —  j  sin  («t*) 


*(•) 


p  i  f  N 

"  ^  ^  —  ^)e  “r  e-^“Zd2  *  ]T^COS  («t») 


ri  W 

Let C(«)  =  [^ cos  (.t*)  — j^sin («t*)] 


where|C(-)j  *  ^ ^2  cos  +  |lsin  («t»)J 


$(«)  “  arc  tan 


(Ssin  (.U)\ 

2  COS  («t»)/ 


The  integral  can  be  evaluated  analytically  and  it  can  be  shown  that 


t. 

JV.  (1  —  e-ie-Klz  -  B(-)e^<“> 


where  B(«)  «  p„ 


t.  M  1  +  (e-t,)2  -2  (ewt.)  sin  (.t.) 


(1  +  <-U2) 


*(•)  »  arc  tan  |  -  2  arc  tan  (»t») 


Thusg(-)  -B(-)  C(-)ej  +  *(“)]  -  |A(-)[ej*(-) 


The  inverse  Fourier  transform  of  p(t)  by  which  p(t)  is  synthesized  from  sinusoids  with 
negative  and  positive  frequencies  is 


P(t)  -  ±  J* gWei^d-  +  *(-)d. 
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The  form  of  the  inverse  transform  which  involves  positive  frequencies  only  is 

2  A(»)  cos(*»t  + 

This  expression  can  be  obtained  from  the  first  form  as  follows . 
p(t)  -  ^  J^KWe^d-  ~  ^ 

P(t)  -i 

p(t)  -  ~  J*  [2(-)ei-t+  «(— )e“j‘*t]  d. 

Where  the  asterisk  denotes  the  complex  conjugate.  Substituting  for  g  (.)  and  g  (-•)  in  the 
last  expression  for  p  (t) 

,<„  -  W  +  ♦«!  ♦  [J  +  #(-)1]o- 

p(t)  -  ^  C 2IA(.)(coe  [.t  4-  *  (.)]  d. 

t 

P(t)  -  J%|B(.)!  |C(«)|  COS  [«t  +  «(-)  +  *  (•)]  <*- 

The  amplitude  spectrum  of  p(t)  is  the  continuous  function  of  frequency  i*  B(.)  C  (-),  .>0. 
B(„)  and  *(•)  can  be  evaluated  from  the  expression  given  above  for  any  specified  value 
of  t*. 

C(.)  and  #<•)  can  be  evaluated  if  the  N  values,  U,  are  known.  They  will  be  determined 
for  the  profile  shown  in  Figure  3. 
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The  twth  respective  times  at  which  the  leading  edge  of  the  pulses  occurs,  are  determined 
from  the  recursion  equation 

U-t^  +  fCUO.k-^S, . .,N 

t.  being  known.  The  value  of  t,  affects  the  phase  spectrum  but  not  the  amplitude  spectrum 
of  p(t).  The  amplitude  spectrum  is  primarily  of  interest  so,  for  convenience,  it  henceforth 
will  be  assumed  thst  t,  =  0. 

f(t*,)  is  a  single  valued  function  of  time  whose  value  is  the  period  of  the  (k-l)th  pulse. 
The  recursion  equation  thus  states  that  the  time  at  which  the  leading  edge  of  a  pulse  occurs 
equals  the  occurrence  time  of  the  leading  edge  of  the  previous  pulse  plus  the  period  of  the 
previous  pulse. 

Figure  3  shows  an  idealized  profile  of  the  period  of  a  pressure  pulse  as  a  function  of  the  time 
of  occurrence  of  its  leading  edge.  It  also  can  be  considered  to  be  the  profile  of  the  firing  rate 
as  a  function  of  the  time  that  has  elapsed  since  the  start  of  firing. 

This  profile  consists  of  three  parts  corresponding  to  the  time  intervals  during  which  the  pe¬ 
riod  is  decreasing  (Part  1),  constant  (Part  2),  and  increasing  (Part  3). 

In  Figure  4,  it  is  assumed  that  K,  L,  and  M  pulses  occur  during  Parts  1, 2,  and  3,  respectively, 
where  K,  L,  and  M  are  integers  and  N  *»  K  +  L  +  M. 

The  recursive  equation  applicable  to  each  part  of  the  profile  is  shown  as  follows: 

Part  1 :  t,  ■  o 

t'  ■  t*-i  t  (2  —t^i)  T,  q  “  2, 3, . . . ,  K 
ts-.-.t.-t-t-  (2-t*)T 


Part  2:  t,  -  V,  +  (2  -  t.)T,q  -  K  +  2,  K  +  3 . K  +  L 

+  (2  —  t,)T 


Part3:  t,  -  V,  +  (2  -  t.)T  +  (V,—  t*)T 

q-K  +  L  +  2,K  +  L  +  3 . K  +  L  +  M-N 


From  these  equations  it  can  be  seen  that  the  period  of  the  first  pulse  in  Part  1,  t,— t,,  is  2T 
and  that  the  period  of  all  pulses  in  Part  2  is  (2  —  t.)  T  which  would  equal  T  for  t,  ■*  1. 
The  period  of  the  first  pulse  in  Part  3  is  (2- 1.)  T  and  the  period  of  the  last  pulse  is  approx¬ 
imately  2T. 

If  the  sequence  of  N  pulses  is  assumed  to  be  periodic  with  a  period  of  T  -  3  see  instead  of 
being  non-periodic  (as  was  conveniently  assumed  in  the  preceding  analysis  and  derivation 
of  the  Fourier  transform,  |g(.)|)  then  the  coefficient  |a*|,  of  the  corresponding  Fourier 
scries  may  be  expressed  in  terms  of  the  transform  as  follows: 
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And  the  magnitude  of  the  spectra  is: 


*w 


The  magnitude,  2jakl,  from  the  foregoing  equations  has  been  determined  with  a  com¬ 
puter  program  for  an  average  firing-rate  of  75  rounds/sec  over  1.0  sec,  a  steady  firing- 
rate  of  100  rounds/sec  for  1.0  sec,  and  a  final  average  firing-rate  of  75  rounds/sec  over 
1.0  sec.  The  line  spectrum  magnitude  of  these  coefficients  as  a  function  of  the  harmonics 
of  the  variable  firing  rate,  is  shown  in  Figure  5.  A  check  of  the  magnitude  of  this  co¬ 
efficient  at  frequencies  other  than  the  harmonics  indicated  the  maximum  value  occurred 
at  each  harmonic  as  shown  in  Figure  6.  Thus,  it  is  seen  that,  for  this  model  of  the  ex¬ 
citation  from  a  variable  firing-rate  gun,  the  Fourier  coefficient  has  a  maximum  value  at 
the  harmonics  of  the  constant  firing-rate  portion  of  the  model. 
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REGRESSION  ANALYSIS 

Regression  analysis  is  the  determination  of  the  position  of  a  line  which  yields  a  minimum 
value  for  the  sum  of  squares  of  the  deviations  of  the  data  values  from  the  line  in  the  y 
(dependent  variable)  direction.  This  line  represents  the  locus  of  the  average  values  of 
the  dependent  variable.  The  position  of  the  line  is  determined  in  terms  of  its  slope  and 
intercept  utilizing  the  method  of  least  squares. 

A  typical  set  of  data  plotted  in  the  Cartesian  plane  is  shown  in  Figure  19.  An  examina¬ 
tion  of  this  plot  indicates  that  the  data  could  be  fitted  by  a  curve  with  the  functional 
form 


G«(f)  -  j8D  «  (28) 


Since  the  solution  of  the  equations  used  in  the  method  of  linear  least  squares  requires 
that  they  be  linear  in  the  unknown  parameters,  it  is  necessary  to  transform  equation  28 
to  a  linear  form  by  taking  the  logarithmic  function  of  each  side  of  the  equation  as  follows 


logG„(f)  =log/J  — zlogD  (29) 


To  simplify,  let 

x  ■*  log  D 

y  -  log  G„  (f) 

•  ■•log/? 

y  "  «-xx  (30) 

Equation  80  represents  the  curve  in  the  log-log  plane  and  is  of  a  form  amenable  to  a  re¬ 
gression  analysis  by  the  method  of  linear  least  squares.  The  sample  of  test  data  shown  in 
the  Cartesian  plane  in  Figure  19  is  spread  out  and  more  easily  analyzed  when  plotted  in 
the  log  G, ,-log  D  plane  as  shown  in  Figure  20.  The  true  values  a  and  z  of  Equation  30  are 
replaced  by  the  estimate  a  and  b,  respectively,  in  the  regression  equation  of  a  line 
through  the  data  points : 

y*-a  +  bx  (81) 

The  properties  band  a  are  defined  as  follows:  (82) 

2  (x*  -  x)  (yk  -  y) 

b-  Vr1 -  (33) 

2  (x*  -  x)* 

k«1 
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where 


a  *“  y  -  bx 

a  m  intercept  of  the  regression  line 
b  m  slope  of  the  regression  line 
x  *  t  -ithometic  avenge  of  the  log  D  data  values 
y  *  arithemetic  average  of  the  log  G0(f)  data  values 

The  computer  program  used  in  this  regression  analysis  was  obtained  from  Reference  11. 

In  this  analysis,  a  confidence  level  (1-a)  of  0.95  was  selected,  corresponding  to  a  signifi¬ 
cance  level  of  0.05. 

After  the  desired  confidence  level  (1-a)  has  been  selected,  it  is  then  possible  to  determine 
the  confidence  interval  ylb  within  which  the  estimated  line  as  a  whole  may  fall.  The 
upper  and  lower  limits  are  determined  for  selected  values  of  x  from  the  following  equa¬ 
tion  given  in  Reference  12. 


y«» 


y  +  b  (x  -  x)  ±  V  2P77 


8/ 


(x  ~  x)f 
“fir  _ 

X  (x*  "*  x)* 

k— 1 


(34) 


where 

F,.«.fdr  (2,  n-2)  degrees  of  freedom  is  found  from  Table  A-5  (T-7)  of  Reference  2. 
F  ”  distribution  of  variance,  joint  statistic 
S,'  —  standard  error  of  the  estimate 


ffi - 

/x(y„~  a-bx*)* 

- 


(35) 


The  index  of  the  goodness  of  fit  or  multiple  correlation  coefficient  is  obtained  from  the 
following  equation: 


lx 

C.-  -  (86) 

A  regression  analysis  was  conducted  on  the  vibration  test  data  from  the  gun  configurations 
SUU-16,  F-4,  MK-4,  and  F-6.  The  results  of  this  analysis  are  listed  in  Table  IV  and  plotted 
in  Figure  22,  23,  25,  26,  27, 87, 88, 89,  and  40.  The  slope  and  intercept  of  the  regression 
line  represented  by  Equation  81  were  obtained  from  Equations  82  and  33.  The  upper  and 
lower  limits  of  the  confidence  band  about  the  regression  line  are  computed  from  Equa¬ 
tions  84  and  85  and  are  plotted  as  slight  curves  in  Figures  22,  23,  25,  26,  27,  37,  88,  39, 
and  40.  The  standard  error  of  the  estimate  and  the  correlation  coefficient  of  each  regres¬ 
sion  analysis  were  obtained  from  Equations  85  and  86  and  are  listed  in  Table  IV. 
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TABLE  IV 


DATA  FROM  REGRESSION  ANALYSIS 


Regression  Line 

Configuration 

Slope 

z 

Standard  Error 
of  Estimate  S,' 

Correlation 
Coefficient  Cf 

UNNORMALIZED 

SUU-16 

-3.70151690 

.7782 

.57883 

F-4 

-3.16325936 

.3398 

.84629 

MK-4  (N  -  42) 

-5.89516664 

.4648 

.88419 

MK-4  (N  *  49) 

-6.39002055 

.4749 

.89391 

F-5 

-1.17115337 

.5058 

.45514 

POOLED  DATA, 
UNNORMALIZED 

SUU-16,  F-4 

MK-4,  F-6 

-2.64185137 

.7062 

.61622 

SUU-16,  F-4, 

MK-4 

-4.67850122 

.6806 

.71425 

POOLED  DATA, 
NORMALIZED 

SUU-16,  F-4, 

MK-4,  F-6 

-2.86865619 

.6881 

.65719 

SUU-16,  F-4, 

MK-4 

-4.55262142 

.6676 

.71900 
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APPENDIX  III 

GUNFIRE  VIBRATION  DATA 

The  data  samples  (Figures  41  through  44)  show  the  truncated  response  function  (de- 
sorilied  in  Section  IV)  superposed  over  the  rumples.  Note  that  best  agreement  results 
when  D  <  76  inches  for  the  F-4,  SUU-16,  and  MK-14  samples.  This  seems  to  be  so  be¬ 
cause  the  gunfire  levels  rapidly  converge  to  the  ambient  vibration  levels  as  D  increases. 
As  this  happens,  it  apitears  that  the  low-frequency  gunfire  contribution  does  not  decay 
as  rapidly  as  the  high  frequency  content,  thus  resulting  in  a  gradually  modified  spec¬ 
tral  shape. 

The  F-r»  data  seems  to  reflect  the  influence  of  the  deflectors,  resulting  in  a  marked 
alteration  of  the  spectrum  from  that  normally  expected.  The  consequences  of  the  de¬ 
flector  influence  are  discussed  in  Section  VI. 
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Figure  41.  F-4  Gunfire  Vibration  Lord 
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Figvr*  42.  SUU-1*  Gunfire  Vibration  Level 
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Figure  44.  F-5  Gunfire  Vibration  Level  (Continued) 


MIL-STD-810B 
NOTICE  2(USAF) 
29  September  1969 


APPENDIX  IV 

MILITARY  STANDARD  ENVIRONMENTAL  TEST  METHODS 


TO  ALL  HOLDERS  OF  MIL^STD-810B: 

1.  The  following:  new  method  is  to  be  added: 

NEW  METHOD  DATE 

T519  Gunfire  Vibration,  Aircraft  29  September  1969 

2.  This  U.S.  Air  Force  Notice  forms  a  part  of  MIL-STD-810B  dated  16  June  1967.  The 
general  requirements  of  MIL-STD-810B  shall  apply  to  any  application  of  this  notice. 

3.  THIS  NOTICE  IS  NONCUMULA'T’TVE,  RETAIN  AND  INSERT  BEFORE  THE  TA¬ 
BLE  OF  CONTENTS. 

4.  Holders  of  MIL-STD-810B  will  verify  above  addition  has  been  entered.  This  issuance, 
together  with  appended  pages  is  a  separate  publication.  Each  notice  is  to  be  retained  by 
stocking  points  until  the  Military  Standard  is  completely  revised  or  cancelled. 


Custodian:  Preparing  activity: 

Air  Force — 11  Air  Force — 11 

(Project  No.  MISC-F650) 


FSCMISC 


MIL-STD-810B 
15  June  1967 


This  tentative  test  method  has  been  approved  by  the  Department  of  the  Air  Force.  It  is 
optional  for  use  by  all  activities. 


METHOD  T519 

GUNFIRE  VIBRATION,  AIRCRAFT 


1.  Purpose.  The  gunfire  vibration  test  is  conducted  to  simulate  the  relatively  brief  but  very 
intense  vibration  fields  resulting  from  blast  pressure  fields  generated  by  repetitive  firing 
guns  mounted  in,  on,  or  near  the  aircraft  structure. 

2.  Apparatus.  Vibration  equipment  with  required  instrumentation. 

3.  Criteria  for  application.  The  most  severe  vibration  field  results  from  blast  pressure  pulses 
coupling  to  the  aircraft  structure  and  inducing  vibration  fields  that  are  of  maximum 
intensity  near  the  gun  muzzle  region.  This  field,  designated  by  a  curve  A  on  figure  519-1,  de¬ 
creases  in  an  inverse  relationship  to  the  distance  as  viewed  in  a  forward  direction  from  the  gun 
muzzle.  The  vibration  field,  designated  by  curve  B  on  figure  519-1,  decays  similar  to  the 
primary  field  beginning  at  the  gun  muzzle  and  extending  in  the  aft  direction.  This  field, 
beyond  the  near  field  region  of  the  muzzle,  is  of  lower  intensity.  In  no  case  should  this  test 
method  be  substituted  for  conventional  vibration  tests.  Guns,  physical  locations,  and  ballistic 
parameters  should  be  carefully  and  accurately  identified  prior  to  application  of  this  test 
method.  If  the  maximum  test  spectrum  level  of  the  gunfire  configuration  is  equal  to  or  les3 
than  .04  g2/Hz,  the  gunfire  method  need  not  be  conducted. 

3.1  Sensitive  equipment.  Equipment  found  most  susceptible  to  gunfire  are  those  equipments 
that  are  usually  located  within  a  3-foot  radius  of  the  gun  muzzle  and  are  mounted  on  the 
structural  surface  exposed  to  the  gun  blast.  Prime  examples  are  UHF  antennas  of  the  blade, 
V,  and  the  flush-mounted  configurations,  including  their  brackeiry,  coaxial  connectors,  and 
lines.  Next  in  order  of  failure  susceptibility  is  equipment  mounted  on  drop-down  doors  and 
access  panels,  equipment  mounted  in  cavities  adjacent  to  the  aircraft  surface,  end  finally, 
equipment  located  in  the  interior  of  the  vehicle  structure.  Typical  vulnerable  equipment  in 
these  latter  categories  are  auxiliary  hydraulic  and  power  units  (including  mounting  brack- 
etry),  switches,  relays,  IR,  photographic,  communication  and  navigation  equipment  and 
radar  systems,  including  items  either  shock  or  hard  mounted. 

3.2  Determination  of  test  levels 

3.2.1  Selection  of  the  maximum  test  level.  The  maximum  gunfire  vibration  level  is  defined 
by  the  following  equation: 

Gm.,  -  E1/E„(G„/r,jx,)  (80/W*)2  (r,nj;  in  PSD  units  (ga/Hz) 
where:  80/W*  —  1,  for  equipment  weight  W*  <  80  pounds 
and:  80/W*  «  Vi  for  equipment  weight  W*  >  160  pounds 

The  test  level  (G».,)  is  obtained  from  the  basic  normalized  level  (G„/r,ju)  which  is  ad¬ 
justed  to  the  muzzle  energy  (E,),  the  firing  rate  (r,),  and  the  number  of  guns  (n,)  of  the 
particular  gun(s)  under  consideration.  .The  normalizing  energy  (E„)  is  equal  to  39,600  foot¬ 
pounds.  The  last  term  (80/W*)  represents  a  mass  adjustment  factor  and  permits  the  reduc¬ 
tion  of  Gm,  when  the  equipment  weight  is  greater  than  80  pounds.  Table  519-1  designates 
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the  muzzle  energies,  the  energy  adjustment  factors,  and  the  nominal  firing  rates  of  com¬ 
monly  used  gun  configurations.  Figure  519-1  provides  G./rjn,  as  a  function  of  the  distance 
parameter  (D).  The  mass  adjustment  factor  is  determined  from  figure  519-2. 

3.2.2  Determination  of  the  distance  parameter  (D).  The  distance  parameter  represents  the 
vector  distance,  measured  (or  estimated)  from  the  gun  muzzle  to  the  mean  distance  between 
equipment  support  points.  Where  equipment  support  points  are  indeterminate,  the  equipment 
center  of  gravity  shall  represent  the  terminal  point  of  D.  The  vector  D  is  usually  generated 
from  the  orthogonal  distances  referenced  to  the  fuselage  station,  the  water,  and  the  butt  line 
data.  The  D  vector  and  the  computation  is  shown  on  figure  519-3. 

3.2.3  Multiple  guns.  For  configurations  involving  multiple  gura,  the  origin  of  D  is  determined 
from  the  centroidal  point  of  the  gun  muzzles.  Figure  519-1  shows  the  origin  location  for  a 
typical  four-gun  staggered  array. 

Example  No.  1 — A  2-pound  UHF  antenna,  is  to  be  flush  mounted  to  the  under-surface 
of  the  aircraft  and  located  forward  of  the  muzzle  of  four  M-61  type  cannons.  The  expected 
firing  rate  is  25  Hz  (r,  -  25) .  (If  the  firing  rate  is  unknown,  the  maximum  specified  rate 
shall  be  used) .  Since,  in  this  configuration,  each  barrel  fires  independently  and  the  number 
of  guns  is  equal  to  four,  then  n,  ■>  4  (see  1/). 

Steps 

a.  The  antenna  location  determined  from  the  aircraft  spatial  coordinates  (see  figure 
519-3)  is  found  to  be  39  inches  or  D  *  89  inches. 

b.  Refer  to  Figure  519-1.  Select  curve  A  and  obtain  the  ordinate  value  of  G./rju  for 
D  —  89  inches.  In  this  case  G./r«n«  “  3.0  x  10* 

c.  From  table  519-1  (  select  E,/E,  for  the  M-61  gun  (see  tj). 

Here,  E,/E«  ”  1 

d.  Since  the  antenna  weight  is  less  than  80  pounds,  the  mass  adjustment  factor  becomes 
unity  and  a  reduction  of  GM  is  not  allowed  (see  figure  519-2) 

e.  Determine  the  normalized  level,  Gi/r,nt,  by  adjusting  G0/r,zu  by  Et/E,  as  follows: 
G./r.n,  »  E,/E,  (G»/rja<,)  orG,/r,n,  “  (1)  3.0x19* 

f.  Multiply  G,/r,n,  by  n,r,  to  obtain  the  maximum  test  level,  G»„: 

G»i  -  G,/r,n,  (r.n,) 

G«„«S.OxlO*(10*) 

Gm>  *  3.0  g*/Hz 

g.  Apply  Gh,  to  figure  519-5  to  obtain  the  test  spectrum. 

Example  No.  t — An  electronic  equipment  weighing  102  pounds  is  located  in  the  aircraft 
nose  section  forward  of  the  gun  muzzle  The  distance  (D)  is  45  inches.  The  firing  rate  of  the 
revolving  barrels  of  the  cannon  totals  100  Hz  with  n,  =  1. 

The  muzzle  energy  (E,)  is  83,000  foot-pounds. 

Steps: 

a.  Repeating  steps  (a)  and  (b)  as  before: 

Ge/rjio  «  1.9  x  10* 

b.  From  note  2,  divide  E,/E«: 

83  x  10739.6  x  101  -  2.1 
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c.  Adjusting  for  E,/E„: 

Gi/r,n,  -  (2.1)  (1.9  x  1(H) 
or  Gi/r,n,  =  4.0  x  1(H 

d.  Obtain  G„.„  as  before : 

Gmi  -  4.0  x  10-1  (102)  (1)  =  4.0  g*/Hz 

e.  Noting  from  85.1  that  W»  >  80  pounds  and  referring  to  figure  519-2,  where  the  de¬ 
rated  value  (G'„,.,)  is  found  to  be  2  dB  down  (-2  dB)  from  4.0  g*/Hz 

From  dB  power  tables : 

-2  dB  -  LOG,,  G',,,.,/4.0 

or  G'M5  =  2.5  g*/Hz 

f.  Apply  G'm„  to  figure  519-5  to  obtain  the  test  spectrum. 

1/  Some  M-61  gun  configurations  (see  MK-IV,  gun  pod,  table  519-1)  feature  two  barrels 
per  pod;  firing  simultaeously.  In  such  a  case,  n,  -  2  per  pod. 

2/  If  E,  is  unknown,  determine  E,  from  ballistic  data  using  E,  ™  m,v,J/2  where: 

m,  «  the  projectile  mass  and  v,  =  muzzle  velocity.  Divide  by  E„  (see  table  519-T'  and  pro¬ 
ceed  as  in  step  (e). 

4.  Test  procedures.  Test  procedures  used  shall  be  as  specified  in  the  equipment  specification 
or  test  plan. 

4.1  Teat  item  operation.  Unless  otherwise  specified,  the  test  item  shall  be  operated  during 
application  of  vibration  so  that  functional  effects  caused  by  these  tests  may  be  evaluated. 
When  a  test  item  performance  test  is  required  during  vibration  and  the  time  required  for 
the  performance  is  greater  than  the  duration  of  the  vibration  test,  the  performance  test  shall 
be  abbreviated  accordingly.  At  the  conclusion  of  the  test,  the  test  item  shall  be  operated  and 
the  results  compared  with  the  data  obtained  in  accordance  with  section  3,  General  Require¬ 
ments,  paragraph  3.2.1.  The  test  item  then  shall  be  inspected  in  accordance  with  section  3, 
General  Requirements,  paragraph  3.2.4. 

45  Mounting  techniques.  In  accordance  with  section  3,  General  Requirements,  paragraphs 
3.25,  the  test  item  shall  be  attached  to  the  vibration  exciter  table  by  its  normal  mounting 
means  or  by  means  of  a  rigid  fixture  capable  of  transmitting  the  vibration  conditions  spec¬ 
ified  herein.  Precautions  shall  be  taken  in  the  establishment  of  mechanical  interfaces  to 
minimize  the  introduction  of  undesirable  responses  in  the  test  setup.  Whenever  possible,  the 
test  load  shall  be  d  atributed  uniformly  on  the  vibration  exciter  table  in  order  to  minimize 
effects  of  unbalanced  loads.  Vibration  amplitudes  and  frequencies  shall  be  measured  by  tech¬ 
niques  that  will  not  significantly  affect  the  test  item  input  control  or  response.  The  input 
control  accelerometer  (8)  shall  be  rigidly  attached  to  the  vibration  table  or  to  the  interme¬ 
diate  structure,  if  used,  at  or  as  near  as  possible  to  the  attachment  point  (s)  of  the  test  item. 

6.  Test  procedures 

5.1  Procedure  I,  Random  vibration  test  The  test  item  shall  be  subjected  to  random  vibration 
along  each  mutually  perpendicular  axis.  Test  times  shall  be  in  accordance  with  time  schedule 
I  from  table  519-11.  The  instantaneous  peaks  of  the  random  vibration  acceleration  may  be 
limited  to  2.5  times  the  rms  acceleration  level  The  power  spectral  density  of  the  test  level 
control  signal  shall  not  deviate  from  the  specified  requirements  by  more  than  +40  —30 
percent  ( ±1.5  dB)  below  500  Hz  and  +100  —50  percent  (±3.0  dB)  between  500  and 
2,000  Hz,  except  that  deviations  as  large  as  +300  —75  percent  (±6  dE)  shall  be  allowed 
over  a  cumulative  bandwidth  of  100  Hz  between  1,000  and  2,000  Hz  for  equipment  items 
whose  weight  is  equal  to  or  less  than  50  pounds.  For  items  weighing  more  than  60  pounds. 
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the  ."tome  deviation  sliali  be  allowed  over  the  extended  range  from  500  to  2,000  Hz  provided 
the  cumulative  bandwidth  deviation  does  not  apply  to  the  swept  bandwidth  of  noise  as  de¬ 
tailed  in  5.3.I.2. 

5.2  Procedure  //.  Simjle  direction  test.  If  the  equipment  item  is  mounted  with  the  base  pe¬ 
ripherally  attached  to  and  in  the  plane  of  the  aircraft  skin,  then  the  test  direction  may  be 
restricted  to  the  direction  normal  to  the  aircraft  akin  (see  figure  519-6).  The  total  test  time 
in  this  case  shall  be  as  stated  in  time  schedule  II  of  table  519-11. 

5.3  Procedure  III,  Composite  or  alternate  test.  For  locations  near  the  gun  muzzle  and  for 
equipment  weighing  in  the  region  of  60  pounds,  the  test  levels  may  exceed  the  force  capabil¬ 
ities  of  all  but  the  largest  shaker  systems.  For  such  cases,  the  following  test  method  may  be 
substituted  provided  the  following  criteria  are  met : 

when:  G„,„,  >  3.0  g*/Hz  >  8.6  x  1(H  (F,/Wt)* 

where:  W*  >  55  pounds  and:  W,  <  1.2  W* 

F,  =  Maximum  RMS  shaker  force  output — pounds 

WT  -  W,+WA+W, 

W*  *»  Equipment  weight — pounds 

WA  *  Armature  weight — pounds 

Wj  *  Test  jig  weight — pounds 

6.3.1  Selection  of  composite  elements.  A  broad  band,  random  vibration  test  level  G»  shall  be 
selected  equal  to  G,„../4  (6  dB  down).  A  random  noise  signal  of  100  Hz  bandwidth  (3  dB 
down  points)  shall  be  superimposed  on  Ga  having  a  G„,  equal  to  the  original  determined 
value.  The  composite  spectrum  is  shown  on  figure  519-7. 

5.3.1. 2  Test  procedure.  The  test  shall  be  conducted  in  accordance  with  6.1,  except  the  nar¬ 
row  band  noise  shall  be  swept  from  300  to  1,000  Hz,  referenced  to  the  narrow  band  center  fre¬ 
quency.  The  sweep  time  shall  be  in  accordance  with  time  schedule  III  in  table  519-11. 

6.  Test  level  overall  rms.  The  overall  rms  test  level  shall  be  no  less  than  the  area  enclosed  by 

the  solid  curve  of  figure  519-5,  -21  percent  (-2  dB).  The  overall  rms  is  defined  by  the  fol¬ 
lowing  equation.  , 

OAR  -  (1163  Gm.,)T 

7.  Summary.  The  following  details  shall  be  specified  in  the  equipment  specification  or  test 
plan: 

a.  Procedure  number  (see  5) 

b.  Pretest  data  required  (section  3,  General  Requirements,  paragraph  3.2.1) 

c.  Nonoperation  of  equipment  during  test,  if  desired  (see  4.1) 

d.  Unique  or  special  test  considerations. 
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